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Abstract. Recent measurements of the deuteron electromagnetic structure functions 
A, B, and T20 extracted from high energy elastic ed scattering, and the cross sections 
and asymmetries extracted from high energy photodisintegration 7 + 6? — » n + p, 
are reviewed and compared to theory. The theoretical calculations range from 
nonrelativistic and relativistic models using the traditional meson and baryon degrees 
of freedom, to effective field theories, to models based on the underlying quark and 
gluon degrees of freedom of QCD, including nonperturbative quark cluster models 
and perturbative QCD. We review what has been learned from these experiments, and 
discuss why elastic ed scattering and photodisintegration seem to require very different 
theoretical approaches, even though they are closely related experimentally. 
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1. Introduction 

The deuteron, the only A = 2 nucleus, provides the simplist microscopic test of the 
conventional nuclear model, a framework in which nuclei and nuclear interactions are 
explained as baryons interacting through the exchange of mesons. With improved 
nucleon-nucleon force models from the 1990s ||T|], and advances in our understanding 
of relativistic bound state techniques, more accurate calculations of deuteron structure 
are possible. 

During the 1990s there have also been revolutionary improvements in our 
experimental knowledge of deuteron electromagnetic structure. The start of experiments 
at the Thomas Jefferson National Accelerator Facility (JLab) has now made available 
continuous high energy beams, with high currents and large polarization, along with new 
detector systems. Several experiments have now significantly extended the energy and 
momentum transfer range of deuteron electromagnetic studies, including A and ^20 foi" 
elastic ed scattering, and photodisintegration cross sections and polarizations. Existing 
experimental proposals promise to continue this trend. Other laboratories have also 
made several important measurements, generally at lower momentum transfer. 

In this context, a review of the deuteron electromagnetic studies, examining the 
current status of the agreement between experiments and theory, is appropriate. We 
attempt to cover our current knowledge of the deuteron electromagnetic structure, 
focussing on the recent JLab results, and prospects for the future. We do not 
consider experiments that use the deuteron as a neutron target, for example, or for 
studies of the (extended) Gerasimov-Drell-Hearn sum rule, deep inelastic scattering, or 
baryon resonance production in nuclei. Our interest is on experiments that probe the 
conventional picture of a nucleus as composed of baryons and mesons, and that probe 
how far models with these effective degrees of freedom can be extended. Table |l] is 
a summary of some of the JLab exeriments that fit this description, and that we will 
review in the sections below. 

The high precision, large momentum transfer measurements may be sensivitive to 
effects not incorporated in the conventional nuclear model. It seems self-evident that 
probes of short distances, well below the size of the nucleon, should require exphcit 
consideration of the quark substructure of the nucleons. Our review suggests that 
evidence for the appearance of these effects seems to depend on the nature of the 
reaction. In elastic scattering, where only the NN chanel is expicitly excited, a successful 
description is obtained using a relativistic description of the A^A^ channel together with 
a minor modification of the short-range structure of the nucleon current (see Sec. |^). In 
photodisintegration by 4 GeV photons, where hundreds of N*N* channels are explicitly 
excited, an efficient explanation seems to require the explicit use of quark degrees of 
freedom (see Sec. ||). 

This review begins with a survey of deuteron wave functions, and then discusses 
the deuteron form factors, threshold electrodisintegration, and high energy deuteron 
photodisintegration. We also call attention to recent reviews by Gargon and Van Orden 
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Table 1. Some JLab deuterium experiments. 



Experiment Reaction / Observables 



Status 



elastic scattering 
91-026 A 
B 

94-018 A 

^20 



paper published Q 
analysis in progress 
paper published 1 3 1 
paper published |4| 



electrodisintegration 
89-028 recoil proton polarization 

94-004 in-plane response functions 

94-102 high momentum structure 

00-103 threshold d(e,e')pn 



analysis in progress 
analysis in progress 
awaiting beam time 
proposal 



photodisintegration 

89-012 cross sections 

89-019 py, Cx', C,> 

96-003 cross sections 

93-017 cross sections 

99-008 cross sections 

00-007 Py, Cx', Cz> 

00-107 Py, Cx', Cz, 



paper published |5| 
paper published |E:| 
paper published |7| 
analysis in progress 
analysis in progress 
awaiting beam time 
awaiting beam time 



g], and by Sick 0. These reviews contain a discussion of the static properties of the 
deuteron and a survey of recent models of the nucleon form factors, two topics we have 
decided to omit from this work. Both also have an extensive discussion of the deuteron 
form factors. 

2. Deuteron Wave Functions 

Calculations of deuteron form factors and photo and electrodisintegration to the A^A^ 
final state require a deuteron wave function, the final state NN scattering amplitude (if 
the transition is inelastic), and the current operator, all of which should be consistently 
determined from the underlying dynamics. Deuteron wave functions used in the 
conventional nuclear model will be reviewed in this section. 

2.1. Nonrelativistic wave functions 

The nonrelativistic A^A^ wave function of the deuteron can be written in terms of two 
scalar wave functions. In coordinate space the full wave function is 

KU^) = E E ^Yim-rriM xIT k^^m - m,\\m) 

= ^Yooir) xlT + ^T. Y2„.-rnM xlT {21m -m, m,|l m) (1) 
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where Y^^n^ are the spherical harmonics normahzed to unity on the unit sphere, zq = u 
and Z2 = w are the reduced 5* and D-state wave functions, and the + distinguishes this 
from other (relativistic) components of the wave function to be described below. The 
spin part of the wave function is 

\ + )a\+)b 



1ms 
Xab 



-^l\+)a\-)b+\-)a\+)b 



ms = + 
= 
m., = - 



(2) 



\~)a \~)b 

Introducing the familar compact notation for matrix operations on each of the two 
nucleon subspaces 1 and 2 

Aaa' l + = A l+)l , (3) 

where A is any 2x2 operator, we can show that 



^00 X, 



Im 
ah 



CTi • (72 X 



Im 



X 



1ms 



X1^2m-ms(r) Xa™' (2 1 m - m, 1 1 m) 



(ScTi • fcTs ■ r - cTi • 0-2) X 



Im 
12 



(4) 



These identities permit us to write the wave function (|I|) in a convenient operator form 
0: 



47r r 



u[r) oi ■ 02 



w{r) 



(3 (Ti ■ r (T2 ■ r - (Ti ■ 0-2) 



X 



Im 



(5) 



In momentum space the deuteron wave function becomes 

1 



^afem(p) 



i2ny 



(Tre 



^ abm (^) 



u{p) a I ■ a 2 



w{p) 



(3 CTi ■ p 0-2 ■ p - CTi • (T2) 



Im 



(6) 



We use the same notation for both coordinate and momentum space wave functions. If 

u{p) = Zo{p) and w{p) = Z2{p), then 



zi{r) 



2 A"" 



7r JO 



TT JO 



rdr Ze{r) je{pr) 
p^dpzi{p) ji{pr) . 



(7) 



Note the appearance of the factors sJ^/ti, a feature of the symmetric definition (^). 
The normalization condition 



implies 



dr 



p^ dp u^ip) + w'^ip) 



(8) 
(9) 
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Figure 1. Reduced coordinate space wave functions for five models discussed in the 
text: AV18 (solid), Paris (long dashed), CD Bonn (short dashed), IIB (short dot- 
dashed), and W16 (long dot-dashed). 



The D-state probability, 

/■oo 

Pd = / drw^(r) (10) 
Jo 

is an interesting measure of the strength of the tensor component of the NN force, even 



though it is a model dependent quantity with no unique measurable value |^ 



The best nonrelativistic wave functions are calculated from the Schrodinger 
equation using a potential adjusted to fit the A^A^ scattering data for lab energies from 
to 350 MeV. The quality of realistic potentials have improved steadily, and now the 
best potentials give fits to the A^A^ data with a x^/d.o.f ^ 1. The Paris potential 



was among the first potentials to be determined from such realistic fits, and it has since 



been replaced by the Argonne V18 potential (denoted by AV18) ||T3[, the Nijmegen 
potentials and most recently by the CD Bonn potentials The 5* and In- 

state wave functions determined from three of these models are shown in Figs. and 
0. These figures also show S and D wave functions from two relativistic models to be 
discussed shortly. In the right panel of the second figure we plot the dimensionless ratios 
u{p)/us{p) and w{p)/ws{p), where the scaling functions, in units of GeV~^/^, are 
, , 16me 
""'^^^ ^ (me + p2)(l+p2/p2) 

16me p'^/pI 
(me + p^){l +p /Po) 

Here e is the deutreron binding energy and m the nucleon mass (we used me = 940 x 2.224 
MeV^), Po = 0-15 GeV^, and pi = 0.1 GeV^. We emphasize that these scaling functions 
have absolutely no theoretical significance and were introduced merely to remove the 
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Figure 2. Momentum space wave functions for five models discussed in the text (see 
the caption to Fig. 1). The wave functions in the right panel have been divided by the 
scaling functions Us{p) and Ws{p). 



most rapid momentum dependence so that the percentage difference between models 
can be more easily read from the ratio graph. We conclude that the five models shown 
are almost identical (i.e. variations of less than 10%) for momenta below about 400 MeV, 
and that they vary by less than a factor of 2 as the momenta reaches 1 GeV (except 
near the zeros). 



2.2. Relativistic wave functions 



The definition of the relativistic deuteron wave function depends in large part on 
the formalism used to treat relativity. In formalisms based on hamiltonian dynamics 
(discussed in Sec. 3.6.6 ) the wave function in the deuteron rest frame can be taken to 
be identical to the nonrelativistic wave function, and no further discussion is necessary 
until the wave function in a moving frame is needed. In formalisms based on the 
Bethe-Salpeter equation the covariant spectator equation 



18 , or on some other 



quasi-potential equation [|19| ^ , the wave functions usually have additional components 
which do not vanish in the rest frame. 



In the relativistic spectator formalism [|T8[, where one of the two bound nucleons is 
off-mass shell, the wave function is a sum of a positive energy component and a negative 



energy component |TG 



P) = {^a6m(p) Ma(-P, &) + ^afem(p) ^a(P, &) } 



(12) 



where p,6) and fo(p,6) are nucleon spinors for the off-shell particle (particle 2 in 
this case) with Dirac index a. The positive energy part has the same structure as the 
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nonrelativistic wave function, with an S and D-state component. The new negative 
energy wave function has the form 



m-m.(p) Xab" (1 1 m - m, m, 1 1 m) -Vs{p)Yi^{p) xlb 




P + Vs{p) {cTi - as) • p 



Xab ' 



(13) 



where Vt{p) and Vs{p) are two additional P-state components of the wave function, and 

xlb = ^[\ + )a\-)b-\-)a\+)b] (14) 

is the nuclear spin wave function. The names of the P-state wave functions follow 
from the fact that vt couples to the spin triplet {S = 1) and Vg to the spin singlet 
{S = 0) wave function. The equivalence of the two forms given in Eq. (|13D follows from 
identities like those given in Eq. (H). 

Two model relativistic 5* and D-state wave functions were shown in Figs. |l| and 



Both models are based on relativistic one boson exchange model developed in Ref. |21 



Model IIB is a revised version of the model of the same name originally described in 
Ref. and Model W16 is one of a family of models with varying amounts of off-shell 
sigma coupling that were introduced in connection the relativistic calculations of the 
triton binding energy described in Ref. [22|. These models are described further in a 



number of conference talks |^ . The relativistic P-state components are small, but can 
make important contributions to the deuteron magnetic form factor. As Figs. ^ and ^ 
show, the large S and D-state components of these relativistic wave functions are very 
close to their nonrelativistic counterparts. 



3. Elastic Electron Deuteron Scattering 

3.1. Deuteron Form Factors and Structure Functions 

Because of the very small value of the electromagnetic fine structure constant {a = 
e'^/Anhc ~ 1/137), elastic electron-deuteron scattering is described to high precision 
by assuming that the electron exchanges a single virtual photon when scattering from 
the deuteron. In this one-photon exchange approximation elastic scattering is fully 
described by three deuteron form factors |2^, |2^, In its most general form, the 



relativistic deuteron current can be written [2S, 25 



2ml 



+ gm{q') [e(r-g)-r^(e-g)] 



(15) 



where the form factors Gi{Q'^ 



1 — 3, are all functions of Q^, the square of the 



four- momentum transferred by the electron, with q = d' — d and = —q^. [In most of 
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the following discussion we will suppress the explicit dependence of the form factors.] 
In practice, Gi and are replaced by a more physical choice of form factors 



c 



Gi + -r]Gc 



GQ = Gi-GM + {l + r])G:i, (16) 

with Tj = Q'^/4:m% At = 0, the form factors Gc, Gm, and Gq give the charge, 
magnetic and quadrupole moments of the deuteron 

Gc{0) = 1 (in units of e) 

Gq(0) = Qd (in units of e/m^) 

Gm{0) = fJ'd (in units of e/2mrf) . (17) 

The form factors can also be related to the helicity amplitudes of the deuteron 
current (where helicity is the projection of the spin in the direction of the particle three- 
momentum). In the Breit frame (where the energy transfer u is zero) the polarizations 
of the incoming (^) and outgoing (^') deuteron are 

(0,±l,-z,0)/V2 A = ± 

-Q/2,0,0,Do)/ma A = 

(0,Tl,-^,0)/v^ A' = ± 

{Q/2,0,0,Do)/md A' = , 

(where the phases for the incoming deuteron follow the conventions of Jacob and Wick 
for particle 2) and the virtual photon polarization is 

(1,0,0,0) A^ = 0. 




(19) 



Hence, denoting the helicity amplitudes by Gy^^ the three independent amplitudes are 



GUQ') = 2D, [Gc + -vGq 

Gy{Q') = G\iQ') = 2Do (Gc - ^vGo) 

G+o(Q') = -GtiQ') = G-_,m = ~G,4Q') = 2DoVvGm, 



(20) 



where 2Do = \/4:ml + . 

The scattering amplitude in the one-photon approximation is 



m(A;',A')7''m(A;,A) (rfV/^M) 



(21) 



where u and u are electron spinors with A;, A {k\ A') the momentum and helicity of the 
incoming (outgoing) electron, respectively. Squaring (^), summing over the final spins, 
and averaging over initial spins give the following result for the unpolarized differential 
cross section 



da 

dn 



da 

dn 



A{Q') + B{Q')t^n\9/2) 



NS 



da 

dn 



(22) 



NS 
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where SIQ"^, 0) is defined by this relation, and 



da 

dn 



a^E' cos2(0/2) 



E' 

CM- 



1 + — sin^ \e 

rrid 



-1 



(23) 



4EHm\e/2) 

is the cross section for scattering from a particle without internal structure (cta/ is the 
Mott cross section), and 6, E, E', and dQ are the electron scattering angle, the incident 
and final electron energies, and the solid angle of the scattered electron, all in the lab 
system. The structure functions A and B depend on the three electromagnetic form 
factors 

2 



A{Q') = G^Q 



2n , A rr^2\ 



Ac{Q') + Aq{Q') + Am{Q') 



(24) 



where the definitions of Ac-, Aq, and Am should be clear from the context. 

While cross section measurements can determine A, 5, and Cm-, separating the 
charge Gc and quadrupole Gq form factors requires polarization measurements. The 
polarization of the outgoing deuteron can be measured in a second, analyzing scattering. 
The cross section for the double scattering process can be written p6| 

da da 



dVtdVto dQdQ^ 



1 + ^hpxAy sin 02 

+ -7^t2oAzz 7^t2 



(25) 



-^'^20^zz - ^'-21^X2 COS 02 + ^^22(^xx " Ay^) COS 202 

where /i = ±1/2 is the polarization of the incoming electron beam, 02 the angle between 
the two scattering planes (defined in the same way as the shown in Fig. |2^), and Ay and 
the Aij are the vector and tensor analyzing powers of the second scattering. Although 
there is a pz component to the vector polarization, the term is omitted from Eq. ([25| ) as 
there is no longitudinal vector analyzing power; without spin precession, this term can 
not be determined. The polarization quantities Pi and t2me (sometimes denoted T2mi, 
but we will reserve capital letters for target asymmetries) are functions of the form 
factors and the electron scattering angle 



Sp^ = - 

Spz = 
— V2S t2o 



4 
3 



1/2 



V{l + V) GMiGc + hGq) tani^ 



{1 + ri){l + ri sin 



3 
1/2 



GIj tan ^9 sec ^6 



t^GcGq + Iv'G 



l + 2{l+ri) tangle 



^2 



V3St2i = 2r] 
-V3St22 = lvG] 



ri + rj^ sin^ \6 



1/2 



GmGq sec \6 



2 ■/ -M (26) 

The same combinations of form factors occur in the tensor polarized target asymmetry 
as in the recoil deuteron tensor polarization. 

Of these quantities, t2o = ^20 has been most extensively measured; it does 
not require a polarized beam or a measurement of the out of plane angle 02. For 
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measurements of A and T20 at forward electron scattering angles, the Gm terms are 
very small, and one may approximate A and T20 by 

-V2Sf,o = lv GcGq + lrfGl (27) 
Introducing y = 2rjGQ/?>Gc gives 

The minimum of T20 ~ T20 = —a/2 is reached for ?/ = 1. The node in the charge 
form factor, Gc = 0, occurs when S ^rj^G'q, and —\/2ST2q |?7^Gq, giving 
T20 = 

This approximation also makes it clear that T20 largely depends on the deuteron 
structure, rather than the nucleon electromagnetic form factors. In the nonrelativistic 
limit (to be discussed shortly), both Gc and Gq are a product of the nucleon isoscalar 
electric form factor multiplied by the body form factor, which is an integral over products 
of the deuteron wave functions weighted by spherical Bessel functions. Hence, in this 
approximation, the nucleon electric form factor cancels in the ratio y oc Gq/Gc-, and 
T20 depends only on the deuteron wave function. 

We note that the relations above between the form factors and the observables are 
model independent, so it is possible to extract form factors from the data and compare 
directly to theoretical calculations. The most complete form factor determination 
appeared recently in Ref. ||30| (see also the analysis in Ref. [H). We will discuss the 
data below in Sec. Ol after we have reviewed the experiments and the theory. 



3.2. Experimental overview 

The initial measurements of elastic ed scattering were by Mclntyre and Hofstadter in 



the mid 1950s [31|. Since then many experiments have run at several laboratories; 



the fits of Ref. |30| include 269 cross sections from 19 references, dating from 1960 to 



the present||[ Polarization experiments are much more difficult. The first results were 

X An important feature of the recent fits of the world data is that the measured cross sections were refit 
rather than using extracted structure functions or form factors. This is necessary since most extractions 
of A (B) used corrections for contributions of B (A) to their cross sections from earher data. In some 
cases alternative definitions (or incorrect formulas) have been given. A minor point is the definition of 
in some cases the recoil factor E' / E is included, while in our definition Eq. (|2^ ) it is not. The 
magnetic form factor Gm can be in units of e/inid (our convention), e/2mp, or dimensionless, with 
magnitude of 1.714, 0.857, or 1.0, respectively, at = 0, and leading to modified coefficients in Eq. (p^). 
Buchanan and Yearian |Q have an alternate definition of Gc, and Gg, with A cx (1 — r/)^(Gp + Gg). 
Benaksas et al and Galster et al include an extra factor 1 + r/ in the magnetic terms in 

A and B, which changes the dependence of the magnetic form factor, though not its value at 
= 0. Ganichot et al and Grossetete et al [36| b oth use a factor of e^, rather than a^, in their 
definitions of the Mott cross section. Cramer et al [ p7| give a dimensionally incorrect formula for their 
(To (= da I (Kl\is!s)^ with explicitly stated energy factors of E' /{EE') = 1/E, as opposed to our factor 
of E'/E^. 
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published in 1984 and there are now only 20 published ^20 data points, and 19 points for 
other polarization observables. The fits of Ref. include a slightly larger data base, 
with 340 points for momentum transfers up to Q of about 1.6 GeV; this misses only a 
handful of the largest momentum transfer SLAG and JLab data points. 

Forward-angle cross section measurements suffice to determine A, both because 
B is small and because of the tan^^/2 dependence. The magnetic form factor Gm is 
determined from large angle measurements of B, since the A contribution vanishes as 6 
^ 180°. With ^ AEE' sm^{e/2) and E' = E{1 + {2E/md}sin^i9/2))-\ we obtain 
the following relations at ^ = 180°: 

q2 = 4:E^(^l + 2E/md^ ^ (29) 
E =(q' + sjQ^{Q^ + Aml) )/4m,. (30) 

One can see that the beam energies needed for high measurements of B are quite 
low, with E = 0.65, 1.02, 1.35, 1.67, and 1.97 GeV corresponding to = 1^ 2, 3, 4, 
and 5 GeV^, respectively. Note that throughout this review we use Q = \fQ^ to avoid 
confusion with the magnitude of the three momentum transfer q, and we use units of 
GeV and GeV^, not fm~^ or fm~^. 

Accurate measurements require that be known accurately since A and B vary 
rapidly with . Energy or angle offsets of a few times 10~^ could lead to being off 
by up to 0.5%. For both A and 5, this leads to offsets that increase with Q^, reaching 
about 2% at g2 = 1 QeV^ and 4% at = g GeV^. 

While cross section measurements can determine A, 5, and G^/, separating the 
charge Gc and quadrupole Gq form factors requires polarization measurements, most 
often ^20- Goincidence detection of the scattered electron and deuteron, which suppresses 
the background and allows experiments to be performed with moderate resolution, is a 
common technique. 



3.2.1. Experimental status of A Several experiments have measured the structure 
function A at small Q. Of particular note are the high precision, 1 - 2% measurements 
from Monterey |3^, Mainz and Saclay ALS The only measurements at 

moderately large Q are from SLAG ElOl [|lll, Bonn [|7| and GEA plus the two 
recent JLab experiments in Halls A and G 0. Data for several experiments are 
shown in Fig. |^ and summarized in Table see Refs. 0, |^0|] or |Q for more extensive 
listings of data. 

Fig. 1^ reveals an unfortunate history of certain measurements not agreeing to within 
the stated uncertainties. For example, at low Q the Monterey and Mainz data overlap 
well, but the overlap of Mainz and Saclay ALS data indicates problems. The four 
largest Q Mainz points used Rosenbluth separations, with A largely determined from 
forward angles of 50°, 60°, 80°, and 90° at 298.9 MeV. Saclay A data were extracted from 
measured cross sections using previous B data. The closest corresponding Saclay points, 
for the same scattering angles at a beam energy of 300 MeV, have cross section about 
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Table 2. Some measurements of A. Symbols are given for data shown in the figures. 

Experiment Q (GeV) symbol # of Year and 

points Reference 



Stanford Mark III 0.48 - 0.88 □ 5 1965 

Orsay 0.34 - 0.48 x 4 1966 

CEA 0.76-1.15 18 1969 

DESY 0.49 - 0.71 O IQ 1971 

Monterey 0.04-0.14 A 9 1973 

SLAG ElOl 0.89 - 2.00 m 8 1975 

Yerevan 0.12-0.19 not shown* 25 1979 

Mainz 0.04 - 0.39 16 1981 

Bonn 0.71 - 1.14 □ 5 1985 

SaclayALS 0.13-0.84 y 43 1990 

JLab Hall A 0.83 - 2.44 ■ 16 1999 

JLab Hall C 0.81 - 1.34 • 6 1999 



|32| 

|33| 

|42| 

|34| 

|38| 

|41| 

|43| 

|39| 

|37| 

|40| 

|2| 

|3| 



*Have larger errors and are consistent with the other data sets. 



7% smaller; the difference is beyond the quoted experimental uncertainties. Significant 
differences such as this are often obscured by semilog plots or not plotting all data sets. 
The body of data, aside from the lowest Q Orsay point, suggests the correctness of the 
Saclay measurments. Theoretical predictions span the range between the two data sets, 
and do not help to determine which is correct. Thus, a new high precision experiment 
in this range appears desirable. 

The agreement between data from CEA, SLAG ElOl, and Bonn near 1 GeV was 
also unsatisfactory. In discussing these measurements, we will compare to the trend of 
the data as determined by the Saclay and JLab measurements. The GEA data have 
large uncertainties, and are systematically low by about la. This experiment measured 
scattered electrons in a shower counter and deuterons in a spectrometer that used a 
quadrupole magnet with a stopper blocking out the central weak field region. In such 
a case it is difficult to determine the solid angle precisely, and this uncertaintly might 
introduce systematic errors into this data. Alternatively, since the spectra were not 
significantly wider than the elastic peak, it has been suggested that over-subtraction 
of background was a problem. However, the background rates were determined to be 
consistent with expected rates from random coincidences and target cell walls. Bonn 
measured coincidence cross sections at large electron scattering angles. Of. ~ 80° - 140°. 
Using forward angle data from SLAG ElOl, GEA, and Orsay, Bonn determined A and 
B. Slightly inconsistent results from the other experiments led to a small uncertainty 
on the Bonn determination of A. Thus, it is only the largest point, for which there 
was only the large angle Bonn data, that has very significant disagreement with other 
determinations of A. Finally, the lowest SLAG point is high. 

The disagreements between the GEA, SLAG ElOl, and Bonn data were part of 
the motivation for two JLab experiments that determined A. Hall A experiment E91- 
026 measured A for from 0.7 to 6.0 GeV^. Hall G experiment E94-018 [| 
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Figure 3. The data for A at low and moderate Q, normalized as explained in 
Section 3.3.2. The data sets are described in Table ^ Note that the right and left 
panels have different vertical scales. All data referred to in Table ^ are shown except 
the highest Q points from Rcfs. and|4l||. 



measured A in the same kinematics as its t2o points, from 0.7 to 1.8 GeV^. The main 
advantages of these experiments over previous work include the continuous beam, large 
luminosities, and modern spectrometers. The Hall A measurements used > 100 /iA 
beams on a 15 cm cryogenic LD2 target, to achieve a luminosity of approximately 5 
X 10^^/cm^/s, and two approximately 6 msr spectrometers. The Hall C measurements 
used the HMS spectrometer along with the deuteron channel built to measure ^20 with 
the recoil polarimeter POLDER. A feature of this system is that the solid angles of the 
two spectrometers were well matched, to within a few percent. 

In the overlap region, the two JLab experiments show better precision than the 
earlier data and generally good agreement; comparisons of theory to data should 
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Table 3. Some measurements of B. Symbols are given for data shown in the figures. 



Experiment 



Q (GeV) symbol # of Year and 

points Reference 



Stanford Mark II 0.10 - 0.13 

Stanford Mark III 0.48 - 0.68 

Orsay 0.20 - 0.28 

Orsay 0.34 - 0.44 

Stanford Mark III 0.44 - 0.63 

Orsay 0.14 - 0.48 

Naval Research Lab 0.11 

Mainz 0.25 - 0.39 

Bonn 0.71 - 1.14 

Saclay ALS 0.51 - 1.04 

SLAG NPAS NE4 1.10 - 1.66 

JLab HaU A 0.7 - 1.4 



□ 



O 

V 





2 
4 
3 
3 
5 
4 
1 
4 
5 

13 
9 
6 



1964 [53| 

1965 [32 1 

1966 [36| 

1966 [33| 

1967 [54| 
1972 [35| 

1980 |55| 

1981 [39| 
1985 [37| 
1985 [52| 

1987 |g, 111 

unpublished 



*These data sets are not shown (B must be inferred from the publication). 



focus on these results, rather than the older data. However, these measurements also 
show a significant disagreement with each other. Uncertainties in each experiment are 
dominated by systematics of approximately 5 - 6%, with statistical precisions near 1%. 
The Hall C data are systematically larger than the Hall A data by just over 2a, slightly 
over 10%, and there appears to be a tendency of the data sets to diverge with increasing 
momentum transfer. This discrepancy will be decreased by a few percent, but not 
eliminated, by a correction to a lower, more accurate, beam energy in Hall C during 
the experiment. It is unclear if the discrepancy can be further resolved. 

An important experimental point is the use in these experiments, and in many 
earlier ones, of ep elastic scattering to calibrate the solid angle acceptance; a fit to the 
world ep cross section data is often used However, recent high precision polarimetry 
results ^ imply that G^^/G\j is significantly smaller than previously believed, 

with G^^/G\j dropping nearly linearly for Q"^ from about 0.5 to 5.6 GeV^. Refitting 
the world cross section data, with the JLab data for the form factor ratio, decreases 
G^; but enhances G\j by about 2% The new fits imply that the ep cross section 
is generally a few percent larger than would have been calculated previously, less than 
the systematic uncertainties of most experiments, and too small to affect comparisons 
of measurements of the ed cross sections and A. The effects on the theoretical deuteron 
form factor predictions will be addressed below. 

In summary, the structure function A is reasonably well determined up to = 6 
GeV^, if one neglects several poorer data points. There remain regions in which there 
are up to about 10% systematic discrepancies between data of different experiments; 
the resolution of these problems is at present unclear. 



3.2.2. Experimental status of B The highest measurements of the B structure 
function come from SLAG NPAS experiment NE4 [pO| , pl|| , which covered the range 
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Table 4. World data for tensor polarization observables. 



Experiment 



Type 



Q (GeV) Observables Symbol # of Year and 

points Reference 



Bates 


polarimeter 


n Q /I 
U.o4, 


U.4U 


J. 

120 


X 



z 


1 nc /I 


Novosibirsk VEPP-2 


atomic beam 


0.17, 


0.23 


?20 




2 


1985 


Novosibirsk VEPP-3 


storage cell 


0.49, 


0.58 


^20 


A 


2 


1990 


Bonn 


polarized target 


0.71 




^20 





1 


1991 


Bates 


polarimeter 


0.75 


- 0.91 


i20)*21,i22 


• 


3 


1991 


Novosibirsk VEPP-3 


storage cell 


0.71 




220 


* 


1 


1994 


NIKHEF 


storage cell 


0.31 




^20 1T22 





1 


1996 


NIKHEF 


storage cell 


0.40 


- 0.55 


T20 


□ 


3 


1999 


JLab Hall C 94-018 


polarimeter 


0.81 


- 1.31 


t2Q-,t2l-,t22 


■ 


6 


2000 


Novosibirsk VEPP-3 


storage cell 


0.63 


- 0.77 


T20 




5 


2001 



* Not shown in the figures. 
** Preliminary data. 



of 1.20 to 2.77 GeV^. These measurements extended the range of previous data from 
Saclay (which went to 1.1 GeV^), and from Bonn |^ (which went to 1.3 GeV^, 
and gave the results for A discussed above). There is good overlap in all but a few of 
the earliest B measurements. Measurements of B were taken as part of E91-026 for 
= 0.7 to 1.4 GeV^, but analysis is not yet final. 

3.2.3. Experimental status of polarization measurements A summary of the 
world data is shown in Table ^ The first polarization measurements were from 
an Argonne/Bates recoil polarimeter experiment and a Novosibirsk VEPP-2 



experiment ^ using a polarized gas jet target. In the gas jet experiment, a 
polarimeter measured the gas polarization after it passed through the interaction region. 

There were three second generation experiments. An Argonne/Novosibirsk VEPP- 
3 measurement pioneered the use of storage cells, increased the internal target 
density about a factor of 15 over the gas jet alone, and pushed out to 0.58 GeV, near 
the minimum in ^20- Because the polarization of the gas varies in the cell, due to wall 
and beam interactions, it was decided to normalize the gas polarization by setting the 
lowest Q datum, at 0.39 GeV, to theory where the uncertainties are small. Such internal 



targets in storage rings are now common. A Bonn polarized target experiment ||60| had 



large uncertainties. At Bates, the AHEAD deuteron polarimeter was used ||6T1 , |6^] to 
determine ^20 in the range just past the minimum of ^20 to just past the node in Gc- A 
continuation of the Novosibirsk experiment had large uncertainties [^], and was never 



published. Note that, to facilitate comparison between different experiments, the data 
are often "corrected" to an electron scattering angle of 70°, but this adjustment and the 
uncertainty it introduces are small. 



Over the past several years, internal target experiments at NIKHEF p4|, 65] and 
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Novosibirsk have improved the precision of the lower data, over a range of Q ~ 
0.3-0.8 GeV. The improvements in Novosibirsk include higher luminosity resulting from 
an improved atomic beam source and a modified beam tune that allows use of a higher 
impedance storage cell. JLab Hall C E94-018 used the recoil polarimeter POLDER 
to measure to the highest Q^, 1.72 GeV^. 

The overlap of the data is good, but apparent systematic shifts can be seen, as the 
NIKHEF and Bates measurements are more negative than the JLab and Novosibirsk 
measurements; note that this is not a difference between polarized targets and recoil 
polarimeters. The issue of determing at what Gc = is affected by this difference. 



The Bates data pT| , |6^ suggest a larger than do the Novosibirsk |^ and JLab 
data. The fits of Ref. |^ do not include the unpublished Novosibirsk data [^, |66|, and 
average between the Bates and JLab points. 

We do not discuss the data for ^21 and ^22- Because of their dependence on Gm, 
they have not been as useful in providing new information as has ^20- To test time 
reversal invariance, one measurement of the induced vector polarization was made |67]. 
The observed result was consistent with zero. 



3.3. Nonrelativistic calculations without interaction currents 

3.3.1. Theory It is instructive to see how the deuteron form factors are related to the 
free nucleon form factors and the deuteron wave function in the nonrelativistic limit. 
Because the deuteron is an isoscalar target, only the isoscalar nucleon form factors 

E — + ^E 

Gm = Gm + (31) 
will contribute to the form factors. In the nonrelativistic theory, without exchange 
currents or (f/c)^ corrections, the deuteron form factors are 

Gc = G%Dc 
Gq = G%Dq 



G 



M 



2mp 



GI,Dm + G%De , (32) 



where the body form factors Dc, Dq, Dm, and De are all functions of Q^. If we choose 
to evaluate these in the Breit frame, defined by 

qo = iy = 0, Do = ^mj + , rf'^ = {Dq, -|q} , rf''^ = {/^o, |q} , (33) 
then the relativistic and nonrelativistic momentum transfers are identical, = q^, and 
the relativistic nucleon form factors can be used without corrections. Note that, in this 
nonrelativistic limit, only the nucleon electric form factors contribute to the deuteron 
charge and quadrupole structure, while both nucleon form factors contribute to the 
deuteron magnetic structure. 

The nonrelativistic formulae for the body form factors D involve overlaps of the 
wave functions, weighted by spherical Bessel functions 

Dc(Q') = / dr u''ir)+w\r) jo (r) 
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V2ri Jo 



dr 



dr w{r) 



wir) 



2u^(r) 



u[r) 



w'^ir] 



32 {r) 



) jo{t) + V2u{r)w{r) + w'^{r) ^2(1") 



drw 



Mr)+j2{T) 



(34) 



where r = gr/2 = Qr/2. At = 0, the body form factors become 



Dc{0) 
Dq{0) 
DMiO) 
De{0) 



dr 



n^(r) + w^(r) 



50 ^0 

dr 



dr w{r) 



u{r) 



w{r) 



2u^ 



drw'^{r) 



w'^(r) 



2-3P, 



D 



(35) 



giving the nomelativistic predictions 

Qd = Dq{o) 

fid =;u, (2-3 Pd) + 1.5 Pd = 1.7596 - 1.1394 Pd, (36) 

with /is = 0.8798 the isoscalar nucleon magnetic moment. The experimental value of 
the deuteron magnetic moment (in these units) is 1.7139, leading to a predicted D-state 
probability of Pd = 4.0%. However this estimate cannot be taken too seriously because 
the magnetic moment is very sensitive to relativistic corrections and interaction currents 
which can easily alter this result significantly. These contributions will be reviewed 
qualitatively later in this review. 

The study of deuteron form factors is complicated by the fact that they are a 
product of the nucleon isoscalar form factors, G*, and the body form factors, D. The 
dependence of the deuteron form factors on older models of the nucleon form factors 
is well discussed in Ref. 0. A year ago the model of Mergell, Meissner and Drechsel 
|68| (referred to as MMD) gave a good fit, and could have been adopted as a standard. 
Figure ^ shows the MMD isoscalar electric and magnetic form factors divided by the 
familiar dipole form factor 



0.71 



(37) 



(with in GeV^). Note that the MMD model does not differ by more than 20% from 
the dipole over the entire range, suggesting that the dipole approximation works very 
well (on the scale of the experimental errors - see below). However, recent measurements 
of the proton charge form factor are producing a surprising result, and at the time this 
review was being completed the picture was begining to change. 

The recent JLab measurements of both the neutron and proton charge form factors 
now suggest that the isoscalar charge form factor may be well approximated by 

1.91r 



Gs JLab 
E 




5.6r) 



+ (1.0 - 0.1262 Q') ^Pb((5' 



(38) 
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,n3 



c 
o 

c 




0.6 



0.2 



(GeV^) 



Figure 4. Ratios of the nucleon isoscalar form factors: ^^{MMD) / Fd (triple-dot- 
daslied line), G^^(MMD)/ {hsFd) (dot-dashed line), /i^G|;(MMD)/G'^,f (MMD) (Solid 
line), G|;JL^'^ 
(long dashed Une) 



G%-'^''^/Fd (dotted line), Fl/G% (MMD) (short dashed hne), and F{/G% (JLab) 



where Fn is the dipole form factor and r = Q'^/Arrip. This JLab model is a sum of the old 
Galster |3i] fit for the nucleon charge form factor (supported by the recent measurements 
||69| ) and a linear approximation to the new JLab G^j/G^ data (from which the 
charge form factor is obtained by assuming that = fipFd). Figure^ shows that this 
form factor differs significantly from the dipole (and also the previously favored MMD 
model), and may have a significant affect on the theoretical interpretation of the data. 
This will be discussed in Sec. ^.8| below. [The Fi/Ge ratios shown in the figure will be 
discussed in Sec. j.6.6 below.] 

Dividing the individual factors Ac, Aq, and Am [introduced in Eq. (P^] by (G|;)^ 
gives reduced quantities that are (except for the weak dependence on the ratio of 
GmsI I^sGes) independent of the choice of nucleon form factor. The contribution of 
these reduced quantities, which we denote by a^, og, and cm, to the total a = A/{G'^^)'^ 
is shown in Fig. ^. The figure shows that the contribution of the magnetic term, qm, 
to the total a is small for < 4 GeV^ (for most of the range it is less than a few 
percent, reaching 10% at ^ 0.5 and also near 4 GeV^). Above of 4 GeV^ it is 
larger, and very model dependent. This justifies the observation that the A structure 
function can be well approximated by A, as stated earlier in Eq. (|27|) . [Note that the 
new JLab data for G^, discussed briefly above, may enhance the magnetic contributions 
to A above 4 GeV^, but will not change these conclusions qualitatively]. 



3.3.2. Comparison to data How well does this simple nonrelativistic theory explain 
the data? The high data for A provide the most stringent test. In Fig. ^ we compare 
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Figure 5. Plots of a{Q^) and the ratio am{Q^)/a{Q^) discussed in the text. Left 
panel shows contributions for models AV18 and IIB-nonrelativistic: a (solid lines), aq 
(long dashed lines), Oc (short dashed lines), a„i (dotted lines). In these plots the ratio 
G\j/G% is fixed at /ig. In each case IIB decreases more rapidly than AV18 at low . 
Right panel: IIB-RIA (solid line), AV18 (long dashed line), IIB-nonrelativistic (short 
dashed line), all with G\JG% — /is, and AVIS with the ratio of G\j/G% given by 
Eq. (H) (dot-dashed hne). 




4 6 8 2 4 

Q' (GeV^) q2 (GeV^) 



Figure 6. The structure function A for the five nonrelativistic models discussed in 
the text, calculated using the MMD nucleon form factors. The models are labeled as 
in Figs. [1] and [2]. The left panel shows data and models divided by a "fit" described 
in the text. See Table || for references to the data. 

the data for A with calculations using the five nonrelativistic wave functions shown in 
Figs. and The calculations use Eq. (|32D with MMD isoscalar nucleon form factors 
and nonrelativistic body form factors given in Eq. (|34D. In the right panel the data and 
models have been divided by the "fit" described in Eq. (^) below. 
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0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

Q (GeV) Q (GeV) 



Figure 7. Deuteron form factors for the five nonrelativistic models compared to data. 
Tire data for Gc and Gq are from analysis of the complete A and ^20 data sets . 
The data for Gj\/ were extracted from the experimental measurements of B, referenced 
in Table 0. MMD nucleon form factors have been used with the nonrelativistic models. 
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1.5 

A(Q2)/fit 

1.0 

0.5 

0.1 1 

(GeV^) 

Figure 8. The structure function A for the five nonrelativistic models discussed in 
the text. The models are labeled as in Figs. [1] and [2]. See Table ^ for references to 
the data. 

It is easy to see that the nonrelativistic models are a factor 4 to 8 smaller than 
the data for Q"^ > 2 GeV^. Furthermore, since the difference between different deuteron 
models is substantially smaller than this discrepancy, it is unlikely that any realistic 
nonrelativistic model can be found that will agree with the data. If the nucleon isoscalar 
charge form factor were larger than the MMD model by a factor of 2 to 3 it might 
explain the data, but this is also unlikely since the variation between nucleon form 
factor models is substantially smaller than this. [If we use the fit Eq. ( |38[ ) to the JLab 
Gep measurements the discrepancy will be even larger.] We are forced to conclude 
that these high measurements cannot be explained by nonrelativistic physics and 
present very strong evidence for the presence of interaction currents, relativistic effects, 
or possibly new physics. 

A detailed comparison of the nonrelativistic models with the three deuteron form 
factors, Gc, Gm, and Gq is given in Fig. |^. The functions used to scale the data and 
theory in the right-hand panels of the figure are 

1.7487 (,,^^1^,^.^ 

««.H^(e-«".0.01e-«'-)(i,g-'(i,J?l)-\ ,3,) 

where is in GeV^ and me = 0.936 x 0.0022246/0.1972. While some of the factors in 
these expressions are theoretically motivated (note the presence of the dipole form of the 
nucleon form factor) we do not attach any theoretical significance to these functional 
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forms; they merely provide a reasonably simple way to scale out the rapid exponential 
decreases from the form factors. Figure |^ shows that the nonrelativistic models do a good 
job of predicting the form factors to a momentum transfer Q — 0.5 GeV, beyond which 
departures from the data and variations of the models make the agreement increasingly 
unsatisfactory.! 

However, careful comparison reveals that there are still (small) discrepancies 
between the data and the nonrelativistic theory, even at low Q^. The data and curves 
from the lower panel of Fig. |^ are shown on an expanded logarithmic scale in Fig. ^ In 
the lowest range from about 0.15 to 0.4 GeV^ the data lie below the nonrelativistic 
theory, and are larger than the nonrelativistic theory only for above 1 GeV^. The 
very low discrepancy seems to be due in part (but not entirely) to the Columb 



distortion corrections that have been used recently to explain the deuteron radius |44 
We will discuss these corrections in Sec. |3]^ below. 

Before we turn to a detailed discussion of the possible explanations for the failure 
of nonrelativistic models to explain the form factors at high Q^, we discuss the low 
momentum transfer results from the perspective of effective theories. 



3.4- Effective field theory 

The recent development of effective field theory provides a powerful method for 
theoretical study of low physics. We will briefly review these results here, and 
return to the discussion of the high results in the next section. 

Effective field theory techniques exploit the fact that the physics at low energies 
E ^ Mq (or large distances A ^ Aq = I/Mq) cannot be sensitive to the details of 
the interactions at very high energies E ^ Mq (or short distances). For example, a 
low energy long wavelength probe may detect the presence of a small scattering center, 
but cannot resolve its structure (much as the far-field of a collection of electric charges 
depends on only one parameter, the total charge). The parameters that depend on the 
short-range physics may be very important, but they cannot be calculated and must be 
determined by a fit to the data. 

Effective field theory works best if the distance scales of the (unknown) short- 
range physics and the (known) long-range physics are clearly separated. Then for 
energies well below the scale of the short-range physics (which we take to be Mq), 
the short-range physics is treated systematically by expanding in powers of E/Mq. In 
applications to the A^A^ system, two scales have been discussed. The so-called "pionless 
theory" chooses Mq ~ m^r, and therefore requires no theory of the ttN interaction. This 
approach can work only at very low energies. The chiral theory chooses Mq ~ ntp 

§ We are not inclined to take the discrepancy between theory and the first Gg point seriously; 
kinematic factors make it difficult to extract this point accurately and it is only one standard deviation 
from the theory. The large Gq and small Gc values for the points at 0.55 GeV and 0.58 GeV result from 
the ^20 data points, from |65[] and |^9| respectively, being about 1 standard deviation more negative than 
calculations and overlapping the negative limit for ^20 of — V2- Note that the tabulated uncertainty of 
Gq at Q = 0.55 GeV in Ref. g should be asymmetric, +0.075/-0.713, (as shown in the Fig. |). 
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and attempts to describe NN scattering up to the p mass scale using the known pion- 
nucleon interaction as given by chiral symmetry. (More precisely, if the magnitude of 
the center of mass relative momentum |p| < Mo/2, the nucleon lab kinetic energy will 
Eiab < M^/2m, which is Eiab < 10 MeV for the pionless theory, and Eiab < 320 MeV 
for chiral perturbation theory.) 

The effective range theory introduced by Bethe is an early version of what 
we now call the pionless effective theory. Weinberg first applied modern chiral 
perturbation theory to A^A^ scattering. He proposed making a chiral expansion of the 
A^A^ potential, and then inserting this potential into the Schrodinger equation. Later 
Kaplan, Savage, and Wise (KSW) ||72|] criticized the consistency of this approach, and 
introduced an alternative organizational scheme, sometimes referred to as Q counting, 
in which the pion interaction is to be included as a perturbative correction (as opposed 
to including it as part of the potential, and counting it to all orders, as proposed by 
Weinberg). KSW applied this method to calculation of the deuteron form factors [^. It 
is now known that the tensor part of the one pion exchange interaction is too strong to be 
treated perturbatively, and recent work has focused on how to include the singular parts 
of one pion exchange in the most effective manner |]74| , [75| . In the following discussion 
we review the recent results from Phillips, Rupak, and Savage (PRS), who give a nice 
account of the calculations of the deuteron form factors in a pionless theory . 

The effective Lagrangian density for a pionless effective theory of the A^A^ 
interaction in any channel (the coupled ^Si — ^ Di for example) is 



{^ip^ tp) ^ (y^ip^ ^p + ip^ V^V^ - 2 VV'^ ■ Wip) + h.c. 



+ ■■■, (40) 

where ip is a. (nonrelativistic) nucleon field operator and Co, C2, and the general 
coefficient C2„ (which fixes the strength of the terms with 2n derivatives) are determined 
from data. The coefficients C2n parameterize the strength and shape of the short range 
interaction. The scattering amplitude predicted by is a sum of bubble diagrams 
which can be regularized using the KSW dimensional regularization scheme with power 
law divergence subtraction |7^. In lowest order (LO) this bubble sum is 

-Co(/x) 47r / 1 \ 



M 



1 + —Co(fi)U + ip) ^ \pcotSo - ip, 
An 

47r / Z, 



, -^]+R{p) (41) 
m \ 7 + «p / 

where p is the magnitude of the nucleon three-momentum in the cm. system, R{p) is 
regular at the pole p = 27, 7 = y/me with e the deuteron binding energy, is related 
to the asympotic normalization of the deuteron wave function, and the dependence of 
Co on the (arbitrary) renormalization point is dictated by the requirement that the 
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Figure 9. Predictions of the cliarge and quadruple deuteron form factors from the 
pionless effective theory, as developed by PRS Q . The data were extracted in Ref . |]30| 
and previously shown in Fig. H. The "fits" are from Eq. (39). 



overall result be independent of ji. The LO result is 

Co(/i) = — , Z, = l, R{p) = . (42) 

In terms of the effective range expansion 

pcot(5o = + |Pd(/ + 7') + Mp^ + 7')' + ■ ■ ■ (43) 

with a the scattering length and pd the effective range, the LO calculation gives 

a = 1/7 , pd = 0. (44) 

Contributions from the next to leading order (NLO) term C2 changes the relations 
in (^2p and (0); in particular, the wave function renormalization constant Zd begins 
to differ from unity and the effective range pd to differ from zero. PRS point out that 
the most stable results are obtained by constraining Co and C2 to give the experimental 
values of the deuteron parameters 7"^ = 4.319 fm and Zd — I = 0.690 instead of 7 and 
Pd = 1.765 fm. This is because the asympotic deuteron wave function is fixed by 7 and 
Zd 



^Ur) = f-§'^ (45) 

and it is the wave function and not the scattering that largely determines the deuteron 
form factors and other electromagnetic observables. 

Using this approach, the LO charge form factor is given entirely by the asymptotic 
wave function (ESf) 



Gg^(Q^) = ^tan-^|-, (46) 
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where g = |q| is the magnitude of the three momentum transferred by the electron, and 
working in the Breit frame (where the differences between relativistic and nonrelativistic 
theory is a minimum) is also equal to \^Q^. PRS show that expansion of the charge 
form factor up to NNLO terms is 



GciQ') = gP{Q'') - {Z, - 1) [1 - G'^'iQ 
-lr%Q'GP{Q') + ... 



(47) 



and because the wave function is correctly normalized, there are no wave function effects 
beyond NLO (the second term). At NNLO effects from the finite nucleon size, r^, 
appear. Similarly, the LO quadrupole form factor obtained by PRS is 

'-nry — \ .-^n -t- \xvy i iM.ri 

47 



'Q 



4g7 - (3g2 + 167^ 



2v^7g3 _ 

with Tjsd the asymptotic S/D-state ratio, and the LO quadrupole moment, Q^^- 
to rjsd/y/^'-f'^ = 0.335 fm^. The expansion of Gq to NNLO is then 



(48) 
equal 



.(0) 

1 
6 



Gq{Q') = G'^>{Q')+miAQd + (Z, - 1) G'-'\Q') - m^Qii 



rlQ'G^S\Q') 



+ 



(49) 



Note that the quadrupole moment at NLO includes a contribution AQd from a four- 
nucleon-one-photon contact term, not determined by A^A^ scattering, and is used to 
fit the experimental value of Qd- PRS suggest that the absence of this piece of short 
distance physics in conventional calculations may explain their underprediction of the 
quadrupole moment. The finite size of the nucleon again comes in at NNLO. 

With parameters largely set by other data, the deuteron charge, quadrupole, and 
magnetic form factors are well predicted up to about Q = 0.2 GeV, as shown in Fig. 
The approach seems to converge well, but beyond NNLO more parameters enter, and 
there is less predictive power. The great strength of the pionless effective theory is that 
strips away complexity, revealing the essential physics required to understand the low 
Q results, and showing (for example) the central importance of the asympotic S-state 
normalization Zd- However, as expected, it clearly does not work for Q much beyond 
0.4 to 0.5 GeV. The theory with pions (sometimes referred to as a "pionful" theory) 
will work to higher ||75| . Removal of divergences from these theories is under active 
study. 

We now return to discussion of the reasons for the failure of nonrelativistic theory 
at high Q"^. 



3.5. Alternative explanations for the failure of nonrelativistic models 



In Sec. p.3.2| we showed that the naive nonrelativistic theory cannot explain the deuteron 
form factor data for Q > 0.5 GeV. In this section we classify the possible explanations 
for this failure, preparing the way for detailed discussions to follow in Secs.^]6] and |3I7| . 

The differences between the data and the nonrelativistic theory can only be 
explained by a combination of the following effects 
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(c) (d) (e) (f) 



Figure 10. Exchange currents that might play a role in meson theories, (a) Large 
/ = 1 TT, p, and A currents that do not contribute to the deuteron form factors, and 
(b) possible / = currents that are identically zero. The currents that do contribute 
to the deuteron form factors are shown in the second row: (c) "pair" currents from 
nucleon Z-graphs; (d) "recoil" corrections; (e) two pion exchange (TPE) currents; and 
(f) the famous pi:^ exchange current. 



• interaction (or meson exchange) currents; 

• relativistic effects; or 

• new (quark) physics. 

The only possibihties excluded from this list are variations in models of the nucleon form 
factors, or model dependence of the deuteron wave functions. In the previous section 
we argued that neither the current uncertainty in our knowledge of the nucleon form 
factors, nor the model dependence of the nonrelativistic deuteron wave functions are 
sufficient to provide an explanation for the discrepencies. 

Possible interaction currents that might account for the discrepency are shown in 



Fig. |I0[ Because the deuteron is an isoscalar system, the familiar large 1=1 exchange 
currents are "filtered" out and only J = exchange currents can contribute to the 
form factors. The / = currents tend to be smaller and of a more subtle origin. The 



nucleon Z-graphs, Fig. p!0|c, and the recoil corrections. Fig. |TO|d, are both of relativistic 
origin. (The recoil graphs will give a large, incorrect answer unless they are renormalized 
IffGl , |77| , [78|.) The two-meson exchange currents should be omitted unless the force also 
contains these forces. The famous p7r7 exchange current is very sensitive to the choice 
of p7r7 form factor, which is hard to estimate and could easily be a placeholder for new 
physics arising from quark degrees of freedom. 

In most calculations based on meson theory, the two pion exchange (TPE) forces 
and currents arising from crossed boxes are excluded, and, except for the p7r7 current 
(which we will regard as new physics), the exchange currents are of relativistic origin. 
Additional relativistic effects arise from boosts of the wave functions, the currents, and 
the potentials, which can be calculated in closed form or expanded in powers of (f/c)^, 
depended on the method used. At low calculations may be done using effective field 
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theories (discussed in Sec. p.4| ) in which a small parameter is identified, and the most 
general (i.e. exact) theory is expanded in a power series in this small parameter. In 
these calculations relativistic effects are automatically included (at least in principle) 
through the power series in (f/c)^. Hence any improvement on nonrelativistic theory 
using nucleon degrees of freedom leads us to relativistic theory. 

Alternatively, one may seek to explain the discrepancy using quark degrees of 
freedom (new physics). When two nucleons overlap, their quarks can intermingle, 
leading to the creation of new NN channels with different quantum numbers (states with 
nucleon isobars, or even, perhaps, so-called "hidden color" states). These models require 
that assumptions be made about the behavior of QCD in the nonperturbative domain, 
and are difficult to construct, motivate, and constrain. At very high momentum transfers 
it may be possible to estimate the interactions using perturbative QCD (pQCD). Very 
little has been done using other approaches firmly based in QCD, such as lattice gauge 
theory or Skrymions (but see Ref. |79|). 

We are thus led to two different alternatives for explanation of the failure of 
nonrelativistic models. In one approach the nucleon (hadron) degrees of freedom are 
retained, and relativistic methods are developed that treat boost and interaction current 
corrections consistently. In another approach, quark degrees of freedom are used to 
describe the short range physics, and techniques for handling a multiquark system in 
a nonperturbative (or perturbative) limit are developed. These two approaches will 
be reviewed in the next two sections. While the discussion appears to be focused on 
the deuteron form factors, it is actually more general, and will be applied later to the 
treatment of deuteron photodisintegration. 

Are these two approaches really different? Superficially, of course, the answer must 
be: Yes! However, QCD tells us that all physical states must be color singlets, and 
a basis of states that describes any color singlet state can be constucted from either 
quarks (and gluons) or hadrons (this would not be true if colored states were physical). 
So at a deeper level it appears that either approach (hadrons or quarks) should work, 
and the best choice is the system that can describe the relevant physics more compactly. 
Further discussion of this issue is clearly beyond the scope of this review. 



3.6. Relativistic calculations using nucleon degrees of freedom 

This long section is divided into six parts as follows: (i) Introduction, (ii) Overview of 
propagator dynamics, (iii) Choice of propagator and kernel, (iv) Examples of propagator 
dynamics, (v) Overview of hamiltonian dynamics, and (vi) Examples of hamiltonian 
dynamics. 



3.6.1. Introduction The inhomogeneous Lorentz group, or the Poincare group, is 
described by 10 generators: three pure rotations, three pure boosts, and four pure 
translations. If we require the interactions to be local and manifestly covariant under 
the Poincare group, we are led to a local relativistic quantum field theory with particle 
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production and annihilation . In this case the Poincare transformations of all matrix 
elements can be shown to depend only on the kinematics (i.e. they depend only on the 
masses and spins of the external particles). The disadvantage is that the number of 
particles is not conserved. If perturbation theory can be used, this approach is very 
successful, but in the nonperturbative regime of strong coupling meson theory it leads 
to an infinite set of coupled equations that cannot be solved in closed form. Numerical, 
nonperturbative solutions of field theory can be obtained in Euclidean space for a few 
special cases Methods that limit the intermediate states to a fixed number of 

particles (two nucleons in this case) are more tractable, and all modern calculations are 



based on the choices depicted in the decision tree shown in Fig. |TT. 

In deciding which method to use, if is first necessary to decide whether or not 
to allow antiparticle, or negative energy nucleons to propagate as part of the virtual 
intermediate state. Since nucleons are heavy and composite, so that their antiparticle 
states are very far from the region of interest, some physicists believe that intermediate 
states should be built only from positive energy nucleons, and that all negative energy 
effects (if any) should be included in the interaction. These methods are referred to 
collectively as hamiltonian dynamics and are represented by the left hand branch in the 
figure. Unfortunately, it turns out that this choice precludes the possibility of retaining 
the properties of locality and manifest covariance enjoyed by field theory. Alternatively, 
in order to keep the locality and manifest covariance of the original field theory, other 
physicists are willing to allow negative energy states into the propagators. These 
methods, represented by the right-hand branch of the figure, are referred to collectively 
as propagator dynamics. Including negative energy states tends to make calculations 
technically more difficult and harder to interpret physically, and those who advocate the 
use of hamiltonian dynamics do not believe the advantages of exact covariance justify 
the work it requires. 

Unfortunately, these two methods are so fundamentally different that many 
physicists do not realize that the limitations of one may not apply to the other. For 
example, for some choices of propagator dynamics all 10 of the generators of the Poincare 
group will depend only on the kinematics, and the Poincare transformations of all 
amplitudes can be done exactly. With hamiltonian dynamics this is not the case; some 
of the 10 generators must depend on the interaction, and transformation of matrix 
elements under these "dynamical" transformations must be calculated. Comparison of 
the two methods is therefore very difficult; the language and issues of each are very 
different and one can be easily misled by the different appearance of the results. We 
cannot discuss these issues in detail in this review, and refer the reader to two recent 



references that survey the subject Here we will give a short review of some 

recent calculations, and explain these differences as we go along. 



3. 6.2. Overview of propagator dynamics Propagator calculations all start from the field 
theory description of two (in this case) interacting particles. While some may prefer to 
express the field theory as a path integral, it is also possible to adopt a more intuitive 
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Figure 11. The relativistic decision tree discussed in the text. 



approach and imagine expanding the path intergral as a sum of Feynman diagrams 
(ignoring issues of convergence for the moment). In order to generate the deuteron 
bound state, which produces as a pole in the scattering matrix, it is necessary to sum 
an infinite class of diagrams, written as 

M (p, p'; P) = V(p, p'; P)+ j ^ V(p, k- P)G{k, P)V{k, p'- P) 
fjnv. r d"- k' 

V{p,k-P)G{k,P)V{k,k';P)G{k',P)Vik',p';P) 



(2 



(27r)" 



+ 



(50) 



where V{p,p';P) is the kernel being iterated, G{k,P) the two body propagator, 
Ai{p,p'] P) the scattering amplitude, and the other quantities are defined below. This 
sum is obtained in closed form by solving the integral equation 

/■ d"'k 

M (p, p'; P) = V(p, p'; P)+j V{p, k; P)G{k, P)M {k, p'; P) . (51) 



If the series (50) is compared to a geometric series 1 + z + z"^ + ■ ■ ■, then the solution to 
the integral equation (0) can be compared to the sum of the geometric series 1/(1 — 2;). 
The geometric series converges only when \z\ < 1, but its unique analytic continuation, 
1/(1 — 2;), is valid for all z. Similiarly, it is assumed that the solution to (|5T|) is valid 
even when the series (|50|) diverges. And just as the geometric series has a pole at z = 1, 



the solution to ([5T| ) will have a pole at P^ 



m 



the square of the deuteron mass. 



The amplitudes V, G, and are all matrices in the NN spin-isospin space, and 
are functions of the four-momenta P = Pi + P2 and p = {pi — P2)/2, with pi and p2 
the momenta of the two particles (labeled in Fig. |T2|). The dimension of the volume 
integration is n, normally either 3+1=4 (3 space + one time dimensions) for the Bethe- 
Salpeter method, or 3+0=3 for the quasipotential methods described below. 



If Eq. (^T]) has a homogenous solution at some external four momentum Pq 



the scattering matrix will have an s channel pole (represented in Fig. |T3D, signifying the 
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Figure 12. The summation of ladder diagrams leads to the covariant scattering 
equation. 



Figure 13. The bound state equation holds near the pole in the scattering amplitude. 

existence of a deuteron bound state. The vertex function for the deuteron bound state 
satisfies the equation 

r (P'k 

r(p; Po) = j V(p, k- P)G{k, Po)r{k; Po) (52) 
with covariant normahzation condition 

The covariant bound state wave function is defined by 

^{p;Po) = Gip;Po)T{p;Po). (54) 

One of the advantages of the propagator approach is that the construction of the 
current operator is comparatively straightforward. It follows (at least in principle) from 
summing all electromagnetic interactions with all the consituents everywhere in the 
ladder sum. For bound states described by the Bethe-Salpeter or Spectator formalisms 
(see the discussion below) there are two diagrams, illustrated in Fig. |l^, that can be 
written 

J^{Pi Po) = eJ — ^ vl/t (ip^ _ k^; P^) J^{k[, k^, ku k2) ^ (iPo - Po) 
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(b) 



Figure 14. The current operator follows from the electromagnetic interaction of all 
the consituents in the ladder sum. For bound states, all of these interactions can be 
collected into the two covariant diagrams shown in panel (a), with the interaction 
current, shown in panel (b), constructed from "fcjjNN and ^(j)(t> couplings (where the 
exchanged meson is denoted by 0). 



+ (2^y J^'^^ik'^P;,) nk[,k'„h,h) nk;Po) . (55) 

In the first term, is the sum of the neutron and proton currents [recall Eq. (|3l|) ] 
and we have chosen particle 1 to interact with the photon (always possible because of 
the antisymmetry of the wave function). The interaction current is I'^, and assumes a 
comparatively simple form if the kernel is a sum of single particle exchanges. This case 
is illustrated in Fig. 

Current conservation, 

g^^'^(P;Po) = 0, (56) 

follows automatically |Q from the bound state equation (p^ if the nucleon and 
interaction currents satisfy the following two-body Ward-Takahashi (WT) identities 

q, J^ik[, k',; k2) = \[l + ra] {G-\k- P,) - G-\k'- P^)} 

q, I^{k[, k',; fci, k2) = V(A;', k; P^) i[l + rg] - i[l + rg] V(A;', k; Pq) . (57) 

Note the appearance of ^ [1 + T3] , the isoscalar charge operator in isospin space. The Jn 
identity is the two-body version of the familiar one-body WT identity 

q,fNiK. h) = {S-\k,) - S-\k[)} , (58) 

with Cat = e (0) for the proton (neutron) and the undressed nucleon propagator 
normalized to S'o"^(A;i) = m— ^1. Note that the constraint on the interaction current 
is not zero (and hence the interaction current is not zero) if the kernel depends on the 
isospin or the total four-momentum Pq. 
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3. 6. 3. Choice of propagator and kernel To fully specify a propagator dynamics, one 
must choose a propagator, G, a kernel, V, and current operators Jj^ and Jj. 

Four different progagators have been used in the study of the deuteron form factors. 
The Bethe-Salpeter (BS) equation uses a fully off shell propagator for two nucleons 



Gbs{p,P) = Si{pi)S2{p2] 





^t{f-p) 








2" 




E^ - f ^ - 

^P V 2 


-Pa) 


2" 



(59) 



where Aj(p)/2m = (mj+ j6)/2m is the (off-shell) positive energy projection operator 
and the right hand expression is the propagator for identical particles in the rest frame 
(with P = {W, 0}, and W used as a shorthand for the four-vector P in the rest frame). 
This choice of relativistic equation was the first to be introduced and is perhaps the best 
known. It retains the full integration over all components of the relative four-momentum 
p, and all of the off-shell degrees of freedom (2 for spin x 2 for "p-spin" , where p = + 
are positive energy u spinor states and p = — are negative energy v spinor states) of 
both of the propagating nucleons, for a total of 4 x 4 = 16 spin degrees of freedom. 
The equation has inelastic cuts arising from the production of the exchanged mesons 
(when energetically possible) and additional singularities when the nucleons are off-shell. 
These can be removed by transforming the equation to Euclidean space. However, the 
BS equation, when used in ladder approximation, does not have the correct one-body 
limit. Numerical comparisons of solutions obtained from the sum of all ladder and 
crossed ladder exchanges with ladder solutions of the BS equation, carried out for scalar 
theories, have shown that the ladder sum is inaccurate and that the one-body limit 
requires inclusion of crossed exchanges ||8^. The BS equation has been solved in ladder 
approximation by Tjon and his collaborators |^, and used to calculate the deuteron 
form factors . The fits to the NN phase shifts originally obtained from these works 
are unsatisfactory by today's standards. 

The Spectator (or Gross) equation (denoted by S) [l^ restricts one of the two 
nucleons to its positive energy mass-shell. If particle one is on-shell, the spectator 
propagator is 

Gsip^P) = 2n 54m' -pI) AM) S2iP2) 

2n5(E,-f- po) AtiPi) A^iW - p,) 



2EpW{2Ep - W) 

where pi = {Ep, p}. This has the effect of fixing the relative energy in terms of the 
relative three-momentum so as to maintain covariance and reduce the four dimensional 
integration to three dimensions [n = 3 in Eq. (p])]. Identical particles are treated 
by properly (anti)symmetrizing the kernel. The restriction of one of the particles to 
its positive energy mass-shell also removes the p = — states of one of the nucleons, 
reducing the number of spin degrees of freedom to 2x4=8. A primary motivation 
and justification for this approach is that it has the correct one-body limit, and 
the three body generalization satisfies the cluster property [Q. The equation also 
has a nice nonrelativistic limit that can be easily interpreted. Numerical studies of 
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scalar field theories [|8J] show that the exact ladder and crossed ladder sum is better 
approximated by the ladder approximation to this equation than it is by the ladder 
approximation to the BS equation. The method can be extended to include gauge 



invariant electromagnetic interactions ||83|. Its principle drawback is that the kernel has 



unphysical singularities which can only be removed by an ad-hoc prescription. Results 
from this method will be reviewed in the next section. 

The internal momentum integration can also be restricted to three dimensions in 
such a way that, for equal mass particles in the rest frame, the relative energy is zero 
and the particles are equally off shell. The Blankenbecler-Sugar-Logunov-Tavkhelidze 
(BSLT) equation [|T^ can be defined so that the two propagating particles are on their 
positive energy mass-shell, reducing the number of spin-p-spin variables to 2x2=4. 
However, this equation does not satisfy the cluster property. The approach of Phillips, 



Wallace, and Mandelsweig 12^, ^7j, which we denote by PWM, also puts the particles 
equally off-shell, but includes all negative energy contributions. Setting n = 3 in 
Eq. (^Tl), the PWM propagator for equal mass particles in the cm. system is 



Gpwuip, P) = j dpo Gpwuip, P)= j dpo { Si{pi)S2ip2) + Gdp, P)} 

f A+(p+)A+(p_) Ar(p+)A2(P-) At{p+)A2{P-) Al{p+)At{P-) ] . . 
^ \ {2Ep -W) {2Ep + W) 2Ep 2Ep j ^ ' 

with 4x4 = 16 spin degrees of freedom. This propagator differs from BSLT primarily by 
the presence of the additional Gc term [which contributes the last two terms in the curly 
brackets involving the "mixed" A+A" projection operators] that includes contributions 
from crossed graphs approximately, and correctly builds in the one body limit. The 
retarded kernel to be used with this propagator, in ladder approximation, is 

Vpwm(p, P) = Gp^y,{p, P) 11 rfpo6^PWM(p, P)Vbs{p, P)Gpwm{p, P)] Gp^^{p, P) (62) 

Perhaps the principle obstacle to implementing this method is that construction of 
current operators is problematic, and manifest Poincare invariance is lost (but Wallace 
P^l has recently shown how to compute boosts for scalar particles exactly). Calculations 
using this method will be described in next section. 

3.6.4- Examples of propagator dynamics We now turn to a description of two examples 
of propagator dynamics. 

Van Orden, Devine, and Gross [VOG]. The Spectator equation has been used to 



successfully describe NN scattering and the deuteron bound state and this work 



uses these results to describe the deuteron form factors ||8^ . The relativistic kernel used 
to describe the A^A^ system consists of the exchange of 6 mesons [n, rj, a, 6 (or oq), p and 
uj]. The model includes a form factor for the off-shell nucleon [^, giving a "dressed" 
single nucleon propagator of the form 

h'^{p) (A^— p^)^ 

^diP) = -( Hp) = 772 2\2 , r 2 2A2 ' (^3) 

m— p (A^ — p^)^ + (m^ — p^)"' 
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where is one of the parameters of the modeL Couphng constants and form factor 
masses (13 parameters in all) are determined by a fit to the data and the deuteron wave 
functions are extracted To insure current conservation, the one-nucleon current 

must satisfy the Ward-Takahashi identity 

{p - prj,{p\p) = Hp') {s^\p) - s,\p')} h{p) , (64) 

and this requires an off-shell modification of the single nucleon current. The solution 
used by VOG is 



3,ip,p) = foip,p){F,iQ')r + 



2m 



+ goip',p)F,{Q' 



m— 



,m- 



2m 



-r 



2m 



(65) 



where ^3(0) = 1 but is otherwise undefined [in the applications described below, 
-^3 = Fd where Fd is the dipole form factor of Eq. (pT])], and 



/o(p',p) 
%{p\p) 



While the on-shell form of the current (| 




(66) 



35]) is fixed by the nucleon form factors, and 
the functions /o and go are fixed by the WT identity (|6^), other aspects of the off-shell 
extrapolation of the current (|65|) are not unique. Using this one-nucleon current, and 
recalling that there are no currents of the type shown in Figs. |iy(a) or (b) [we postpone 
discussion of the pTij current], it was shown [Q that the full two body current to 
use with the Spectator equation is given by the diagrams shown in Fig. |15|. These 
diagrams are manifestly covariant, and automatically include effects from Z- graphs or 
retardation illustrated in Figs. p!0|(c) or (d). They are referred to as the complete impulse 
approximation (CIA) to distinguish them from the relativistic impulse approximation 
(RIA), an approximate current used in earlier calculations ||2^. The RIA is obtained 



by multiplying diagram |T5|(a) by two, and is very close to the CIA. 

Phillips, Wallace, Divine, and Mandelsweig [PWM]. This work is based on the 
Mandelsweig and Wallace equation [0, supplemented by contributions from the crossed 
graphs [^, as described above. It is sometimes referred to as the equal-time approach. 
A feature of this equation is that it includes the full strength of the Z-graphs; the 
A+A^ contributions shown in Eq. ( pT] ) are roughly twice as strong as the Z-graph 
contributions included in the Spectator equation. The PWM propagator is also explicitly 
symmetric, a convenience when applied to identical particles. In the published work 
reviewed here [^, the deuteron is described by a one boson exchange force using the 
parameters of the Bonn-B potential with the exception of the a meson coupling, which 
is adjusted to give the correct deuteron binding energy. Lorentz invariance is broken 
by the approximation; the boosts of the deuteron wave functions from their rest frames 
are treated approximately. The PWM current is a modification of Eq. (|55D . In the 
present work retardation effects [like those illustrated in Fig. p!0|(d)] are omitted from 
the current operator; only one body terms and Z-graph contributions are included. 
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(a) 



(b) 



Figure 15. The complete impulse approxiation (CIA) to a ladder Spectator theory. 
The X denotes the on-shell particle. In diagrams (b) the kernel connects the region 
where one of the two particles is on-shell to the region were both particles are off- 
shell. For identical particles, it can be shown that (b) ~ (a), so that the sum of these 
three diagrams is approximately equal to 2x (a), which is referred to as the relativistic 
impulse approximation (RIA). 



We now turn to a discussion of the other major approach detailed in Fig. 0: 
hamiltonian dynamics. 

3. 6. 5. Overview of hamiltonian dynamics Approaches based on hamiltonian dynamics 
start from a very different point than propagator dynamics, and this is one reason 
it is difficult to compare the two. While propagator dynamics starts from field theory 
(which can be described as a quantum mechanics with an arbitrary number of particles), 
hamiltonian dynamics starts from quantum mechanics with a fixed number of particles. 



For a detailed review, see Ref. [B2 



Quantum mechanics begins with a Hilbert space of states defined on a fixed space- 
like surface in four-dimensional space-time. The various options for choosing this space- 
like surface were classified by Dirac in 1949 The instant-form corresponds to 
choosing to construct states at a fixed time to = 0, and is the choice usually made in 
elementary treatments. Alternatively, front-form quantum mechanics constructs states 
on a fixed-light front, customarily defined to be t+ = t -|- z = 0. (We use units in 
which the speed of light, c, is unity.) More generally, the light-front may be chosen 
in any direction defined by x^n^ = 0, with = {l,n} and = 1. Finally, point- 
form quantum mechanics constructs states on a forward hyperboloid, with — = a^; 
if: > [the limiting cases of a = oo gives the instant-form (with to = a), and a = 
the front-form]. These three surfaces are shown pictorially in Fig. |T6|. [While the point 
t = 0, r = is not on the hyperbolid, all distances between points on the hyperbolid 
are space-like.] 

The Poincare transformations are symmetries that leave all probabilities unchanged; 
they must be unitary transformations (with hermitian generators) on the space of 
quantum states. The 10 generators of the full Poincare group are the hamiltonian, 
H, generator of time translations, three-momenta, P*, generators of spatial translations, 
angular momenta, J*, generators of rotations, and K^, generators of boosts. They satisfy 
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Figure 16. Drawing of the three surfaces on which states can be defined in quantum 
mechanics. All points on the forward hyperbolid are separated by a space-like interval. 

the following commutation relations: 

[H, P'] = [H, r] = {P\ P^] = , [J\ X^] = ie,,fcX'= , for X' = J\ P\ K' 

[K\K^] = -ieijkJ', [K\P^] = -i5ijH, [K\ H] = -iP' (67) 

For each of the forms of quantum mechanics there is a subgroup of the Poincare 
transformations that leave the states invariant on the fixed surface associated with that 
form. This is the kinematic subgroup, and the transformations in this subgroup will 
not depend on the dynamics (since the dynamics describe how the states change away 
for the fixed surface). In the instant- form, space translations and rotations clearly leave 
the surface t — invariant. Generators of these transformations form a subgroup of the 
Poincare group, with commutation relations 

[P\ P^] = , [J^ P] = ieijkJ" , [J\ P'] = ieijkP" . (68) 

Transformations of states under these transformations will not depend on the dynamics. 
The hamiltonian H carries the states away from the initial fixed t — surface, and 
contains the dynamics. The other three generators (the boosts) will also, in general, 
depend on the dynamics because their commutators involve H. 

The front-form surface t+ = is left invariant by translations in the x, y, 
and t~ directions [the generator of translations in the t^ direction is because 
Ht — PzZ = {H~^t~ + H~t~^)/2]. It is also left invariant by rotations and boosts in 
the z direction, and by the generalized boosts = + and E^ = — J^. These 
7 generators form a subgroup of the Poincare group, with commutation relations 

[P\P^] = [P\H+] = [r,H+] = [E\E^] = [E\H+] = [r,K'] = [K',P'] = 
[J^ P'] = teijP^ , [J^ E'] = tSijE^ , [E\ P^] = -tSijH+ 
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[K', H+] = -iH+ , [K', E'] = -iE' , (69) 

where i and j = {x,y} = {1,2}, and eu = 1, £21 = —1, and en = £22 = 0. The 
fact that the front-form kinematic subgroup includes seven generators, including the 
boost and generalized boosts E\ makes the front-form popular. But a principle 
motivation for using the front-form is that it is a natural choice at very high momentum, 
where the interactions single out a preferred direction (the beam direction) and the 
dynamics evolves along the light-front in that direction. The disadvantage is that the 
generators that contain dynamical quantities are H~ and J*, and this means that angular 
momentum conservation must be treated as a dynamical constraint. 

Finally, the point-form hyperbolid is left invariant by the homogeneous Lorentz 
group itself, with commutation relations 

[J\ J^] = ie,,kJ^ , [J\ K^] = lei.kK' , [K\ K^] = -le^.^j'^ . (70) 

The hamiltonian and the momentum operators all carry the dynamical information. 

We see that each of the forms of quantum mechanics has a different set of kinematic 
generators, and in no case are they all kinematic. Practitioners of hamiltonian dynamics 
sometimes speak as if it were impossible to treat the full Poincare group kinematically. 
This is true only in the context of hamiltonian dynamics; all of the generators are 
kinematic in the BS or Spectator forms of progagator dynamics. 

Dynamics is introduced whenever the states are propagated away from the surface 
on which they are initially defined. As in normal quantum mechanics, the deuteron will 
be an eigenstate that propagates in "time" without loss of probability; it will be an 
eigenstate of the generalized hamiltonian. In the instant-form, the rest state \E'i(0) is an 
eigenstate of the momentum operators 

P^i(0) = O, (71) 

and the bound state equation in the rest frame is 

i/^i(0) = m,^i(0) . (72) 

In the front-form the rest state \£'f(0) is an eigenstate of the operators P_l = {P^,P^} 
and i/+ 

P^^f(0) = O, /f+^F(0) = mrf^F(0), (73) 

and the dynamical bound state equation is 

if-^F(O) = mrf*F(0) . (74) 

Finally, in point-form the rest frame eigenfunction must satisfy the four dynamical 
equations 

P^^p(0) = P<r^p(0), (75) 

Where = {i/, P*} and = {md,0}. In applications, the dynamical equations 
(0), (|7^ and the /x = component of (|7|) can all be taken to be the nonrelativistic 
Schrodinger equation in the rest frame, so the same nonrelativistic phenomenology can 
be used for any of these forms of mechanics. 
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To complete the calculation of the deuteron form factors using hamiltonian 
dynamics one must choose a current operator that conserves current, and construct 
the proper matrix elements of this operator between deuteron wave functions. This will 
be discussed next. 



3. 6. 6. Examples of hamiltonian dynamics The steps taken to construct the current 
and calculate the form factors depend on the form of quantum mechanics used, and the 
taste of the investigator involved. Here we briefly describe recent work by five groups. 

Forest, Schiavilla, and Riska [FSR]: Based on the work of Schiavilla and Riska ||91|| , 
Forest and Schiavilla have done an instant-form calculation of the deuteron form 
factors. The original work of Ref. |^ used on-shell matrix elements of the one body 
charge and current operators 



E{P')E{P)) 



1/2 



2m 



u{p, s] 



(76) 



expanded in powers of f /c. Here Fi and F2 are the Dirac and Pauli form factors, usually 
replaced by the familiar charge and magnetic form factors 



FiiQ') + F,{Q') . (77) 

In the recent unpublished work of Ref. the calculations have been done in momentum 
space, where the one body current operators have been evaluated without making any 
(w/c) expansions, the relativistic kinetic energy y/m'^ + has been used in place of the 
usual nonrelativistic expansion m + p^/ (2m) [with the parameters of the AV18 potential 
refitted], and the boost corrections to the deuteron wave functions have been included. 
This work also includes two-body charge operators from vr and p exchange using methods 
developed by Riska and collaborators |9^ . 

Arenhdvel, Ritz, and Wilbois [ARW]: This recent calculation ||9^ does a systematic 
v/c expansion of relativistic effects that arise from the one body current operator and 
from contributions from meson exchanges. The current operator (^6]) is approximated 
by 



fAi 



Fi 



(p' + P) 



+ G 



M 



'Am? 
i [a X q 



+ (2Gm - Fi 



la 



iq X P) 



4m^ 



+ 



[v/cf 



/i = 



(7^ 



2m 2m 

where = or (with appropriate Fi and F2), a is the operator in the nuclear spin 
space, and q = p' — p is the t/iree-momentum transferred by the electron. This charge 
operator is correct to order (i;/c)^, and the (f/c)^ contribution to the current operator 
is given in Ref. |^. The q^/8m^ correction term to the charge operator is referred to 
as the Darwin- Foldy term. The o" ■ (q x p) is the spin-orbit term. 

There are ambiguities in all calculations based on expansions in powers of (f/c)^. 
One ambiguity arises from the fact that the square of the three- momentum, q^, depends 
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on the frame in which it is evaluated. In the Breit frame, = Q^, while in the center 
of mass of the final deuteron (the frame preferred by ARW), = (1 + ri)Q'^ (recall 
that rj was defined in Sec. |3]l|). A second ambiguity surrounds the choice of G^; versus 



Fi. Some experts |9T[] advocate using Ge (because it is the correct charge operator) in 
place of Fi. The difference between Fi and G^; is of higher order. These ambiguities 
introduce theoretical uncertainty into any calculation. The size of this uncertainty 
depends on both the value of and the choice of nucleon form factors; for example 
the difference between using Fi or Ge can be inferred from the ratios shown in Fig. ^ 
and is large for the recently measured JLab form factors and small for the MMD form 
factors. Uncertainties of this kind do not arise if the calculation is done covariantly, or 
to all orders in (f/c)^. 

ARW also include boost corrections originally derived to lowest order in v/c by 
Krajcik and Foldy The boost corrections can be written as an operation on the 



wave function of the form 

B^o{r) = e-*^ ^o(r) ~ (1 - ix) ^o(r) , (79) 
where X = Xo + Xv with 

[r ■ P)(p-P) ^ ^ \ ^ [{a, - a^) x p] 



the boost associated with the kinetic energy and the spin and xv the boost associated 
with the potential. In (PP]), r and p are the relative coordinate and relative momentum 
of the nucleon pair, and P is the three-momentum of the moving deuteron. ARW use 
the values of xv worked out by Friar and also include relativistic effects from 



retardation, isobar currents, and meson exchange. To evaluate the latter a meson 
exchange model is needed, and ARW use the interactions and parameters of the Bonn 
OBEPQ potentials [only results from the OBEPQ B potential are presented in the next 
section, although Ref . |Q includes results from all three OBEPQ potentials] . Friar has 
emphasized that relativistic effects can be moved in and out of the wave functions and 



currents by unitary transformations [p7| , |98| , so that all of these effects are ambiguous 



unless fully defined by the theory. Effects due to pair currents or recoil corrections, 
shown in Fig. p!0|(c) and (d), do not appear to be included. ARW state that their 
calculations should be good only up to — 1-2 GeV^. 

Corrections to the charge operator to order (f/c)^ obtained from instant-form 
dynamics and from the Spectator form of propagator dynamics have been compared 
|76| , |97| , |98| , [99|] . In the cases studied, the same total result was obtained from the sum 
of all of the corrections, but the individual terms in the sum were found to have a very 
different form even when they appeared to come from the same physical effects. 

Carhonell and Karmanov [CK]: In this front-form calculation [|100|| the direction 
of the light-front [denoted by cu'^ = {1, n} where = 1] is treated as an unphysical 
degree of freedom. Wave functions and amplitudes may depend on oj but only those 
components of scattering matrix elements independent of uj will be physical. It is argued 
that this approach will give an explicitly covariant front-form mechanics ||101|| . When 
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applied to the deuteron form factors there are 11 spin invariants, three that are physical 
and 8 that depend on uj and are unphysical. In an exact calculation the 8 unphysical 
invariants would be zero, but in approximate calculations, such as that carried out in 
Ref. ||100|| , they will not be zero. The deuteron form factors can be extracted from 



the three physical invariants by projecting them from the general result, as derived in 
Ref. ||101|| . For the choice = {1,0,0, —1} (corresponding to choosing the front-form 



surface t~^ = 0) this method shows that the charge and quadrupole form factors can be 
extracted from the J"*" component of the current (in common with other treatments), but 
also shows that the magnetic form factor cannot be obtained only from this component 
and requires a different projection (and includes contributions from contact terms). The 
rules for a general graph technique for calculating amplitudes in this formalism are given 



in Ref. 102 



Using this method the deuteron wave function will in general have 6 components, 
only three of which have been found to be numerically large. In addition to the familiar 
5" and D-state components, the third large component is proportional to a new scalar 
function /s, and adds the term 

3 i „ , . ,^ ^ Im 




Kbmip) = -xhr^h (^1 - f^2) • (P X n) X 



to the deuteron wave function displayed in Eq. (^) (to obtain this form we renormalized 
the expression in Ref. |p.00|| so that fi = u and used the transformations in Ref. 



In Ref. ||100| /s is calculated perturbatively using the Bonn potential from Ref. | ]103| 



without change of parameters. They find that /s is the largest of the three components 
for all momenta greater than 500 MeV, and believe that the perturbative estimate is 
accurate to about 20%. The physical meaning of the /s contribution has been studied in 
threshold deuteron electrodisintegration, where it contributes about 50% of important 
pair term contributions. 



Lev, Pace, and Salme [LPS]: The LPS ||104|| calculation is a recent version of a series 



of light-front calculations that have assumed the light-front is fixed {dX t^ = t + z = Q). 
In the past, calculations with fixed light fronts have run into a problem with the loss of 
angular momentum conservation, and before we review the LPS results we will discuss 
this issue. 

In calculating form factors with fixed light fronts it has been conventional to choose 
a coordinate system where = g° ± = and g_L = Q. Current conservation is then 
satisfied if only one component of the current {J^) is non zero. Consider the matrix 
elements of the deuteron current, J'^iy^, where A' (A) are the hehcities of the outgoing 
(incoming) deuterons. One consequence of the loss of manifest rotational invariance 
is that there are four independent matrix elements of the current related by the 
constraint 

Jo+ + 2v^J+o - JX- - (1 + 2r/) J++ = . (82) 
This is a dynamical constraint often referred to as the "angular condition" [p.05| , |106|| . 



The deuteron form factors can be extracted from any choice of three of the matrix 
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elements Jy^, and if condition (|^) is not satisfied each choice will yield different 
results. The form factors will not be uniquely determined unless the angular condition 
is satisfied. 

To avoid (or solve) this problem, LPS work in the Breit frame, where = ±Q 
and q± = 0. Current conservation then requires that = J~ . A current operator 
that satisfies these conditions was constructed in Ref. ||107|| . For elastic scattering this 
operator has the form 

1 c I ■ c /■ I \* ■ a ^ \ a = a for a = +, ± , ^ 

J'' = l,{jfree + Lte-'^{J^^^:) e-^-'^} \ . , (83) 

/ V J / J I a = + for a = - 

where J free is the free (one body) current operator, and and e^'^^'^ are rotations by 
— TT about the x axis, in the vector space and e™^"^ in the spinor space. Note that the 
definition insures that = as required by current conservation. Using this current, 
LPS have calculated the deuteron quadrupole moment, Qa, to 2% accuracy ||108|| . The 
calculation shown below in Sec. uses MMD nucleon form factors and the Nijmegen 
II deuteron wave functions. 

Allen, Klink, and Polyzou [AKP]: The deuteron form factors have also been recently 
calculated using the point-form of quantum mechanics | 109 |. Here there is no difficulity 
in writing down manifestly covariant matrix elements, but there is some ambiguity in 
deciding how to impose current conservation. AKP work in the Breit frame, choose a 
one body impulse current to describe the /i = 0, 1, and 2 components of the current, 
and introduce a two body current J| (which need not be calculated) to insure current 
conservation. New effects come from the way the wave functions are constructed in 
point form ("velocity" states are constructed), and from the fact that momentum is 
now a dynamical generator, so that the momentum transferred to the nucleon inside 
the deuteron is not equal to the momentum transferred to the deuteron as a whole. 
They argue that the momentum transferred to each nucleon inside the deuteron is 

Ql = ip[ - Pif = Q'^^^^^^ {i + v)> Q\i + V) . (84) 

At momentum transfers — 4 GeV^ this is a 25% increase, and leads to a large 
suppression of the form factors. This explains part of the decrease in the size of the 
form factors predicted by this model. The results reported below use the MMD nucleon 
form factors and the AV18 NN potentials. 

We now turn to a brief review of methods using quark degrees of freedom. 



3.7. Calculations using quark degrees of freedom 

Calculations based on quark degrees of freedom must confront the fact that the deuteron 
is at least a six quark system. Since the six quarks are identical (because of internal 
symmetries) the system must be antisymmetrized, and it is not clear that the nucleon 
should retain its identity when in the presence of another nucleon. How does the 
clustering of the six quarks into the two three-quark nucleons appear at large distance 
scales? How do we treat the confining forces in the presence of so many quarks? 
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Figure 17. Diagram showing a photon couphng to a quark exchange diagram. 
The exchanged gluon (or meson), shown by the heavy dotted Hne, distributes the 
momentum equally between the two three-quark clusters, suggesting the result ds^). 



The approach to these issues depends on whether of not is large enough to 
justify the use of perturbative QCD (pQCD). 

3.7.1. Nonperturbative methods At modest the momentum transferred by the 
gluons is small and the QCD coupling is too large for perturbative methods to be useful. 
In the nonperturbative regime calculations must be based on models. Many papers have 
been written addressing these issues, and a complete review is beyond the scope of our 
discussion. Here we mention only two contributions that give the flavor of the discussion. 
Maltman and Isgur ||110| , |111|| studied the ground state of six quarks interacting through 
a qq potential previously used to explain the spectrum of excited nucleons, and found 
that there was a natural tendency for the quarks to cluster into two groups of three (i.e. 
nucleons). They obtained a reasonable description of the deuteron, and confirmed that 
the short range A^A^ repulsion could be largely understood in terms of quark exchange. 



Later, de Forest and Mulders ||112|] , using a very simple model, considered the effect 
of antisymmetrization on the structure of the form factor. Their calculations suggest 
that the zeros seen in form factors could be a consequence of antisymmetration alone. 
They also show that the factorized form of the impulse approximation obtained from 
nonrelativistic (and some relativistic) theories, which gives the deuteron form factors as 
a product of a nucleon form factor and a nuclear (or body) form factor 

Gd{Q')=Fn{Q')xD,{Q^), (85) 

may not be a good description in the presence of antisymmetration. When quarks are 
exchanged between nucleons it is no longer possible to separate the nucleon structure 
from the nuclear structure. Consideration of the quark exchange diagram shown in 



Fig. suggests a factorization formula of the form 

= \Fr,{[Q/2Y)] ' x D.iQ') . (86) 



Because nucleons are composite and identical, either of the forms ( p5D or (|86D (or yet 
other relation) might hold, and there is no clearly correct way to isolate the structure 
of the nucleon from the structure of the bound state. 
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In model calculations these issues can be handled by separating the problem into 
two regions: at large separations {R > Re) it is assumed that the system separates into 
two nucleons interacting through one pion exchange, and at small distances {R < Rc) 
the system is assumed to coalesce into a six-quark bag with all the quarks treated on an 
equal footing. In this review we report the results of a calculation by Dijk and Bakker 
|113|| , where references to other calculations of this type can also be found (see also the 



work of Buchmann, Yamauchi, and Faessler ||114|| ) 



Dijk and Bakker [DB]: This calculation is based on the quark compound bag model 
introduced by Simonov | |115| ]. Here the six-quark wave function is assumed to be the 
sum of a hadronic part and a quark part. The hadronic part is a fully antisymmetrized 
product of two three-quark wave functions, each with the quantum numbers of a nucleon, 
and a relative NN wave function Xnn{'^) 

I'iph) = A [^^(123)V'^(456) XNN{r)] , (87) 

where r = (ri + r2 + — r4 — rs — r6)/3 is the effective internucleon separation. The 
quark part is a sum of eigenstates \ipu) of a confined 6-quark system 

IV',) =E«-I^-) > (88) 

V 

where the confined 6-quark states are zero outside of a confining radius r = b, which 
is a parameter of the calculation. In the applications, only one term i> = 1 needs to 
be included in the sum (|88|). The dynamical quantities determined by the calculation 
are the NN wave function XAfAf(r) and the spectroscopic coefficient ai{E) which is a 
function of the energy E. 

The A^A^ scattering phase shifts and mixing parameters are determined by replacing 
the spectroscopic coefficients by boundary conditions on the surface r = b and 
integrating the Schrodinger equation for r > b. The Paris potential [1^ is used to 
describe the A^A^ interaction in the peripheral region and is set to zero in the inner 
region. Two models were developed; in this review we report results from the fits to the 
Arndt single energy 1986 |[L16|| solutions, which DB denote QBC86. This fit finds 6 = 1 
fm. Calculation of the form factors requires an assumption about the form factor of the 
internal compound bag part of the wave function. They use 

F,i(Q2) = (l + 2gV5A?)"' (89) 

with Ai = 1 GeV obtained from a fit to the A and B structure functions. Results from 
this model are reported in Sec. below. 



3.7.2. Perturbative QCD If one believes the momentum transfer is high enough, 
perturbative QCD (pQCD) may be used to study the deuteron form factor and reactions. 
Here it is assumed that the problem naturally factors into a hard scattering process 
in which the momentum transfer is distributed more or less equally to all of the six 
quarks, preceeded and followed by soft, nonperturbative scattering that sets the scale 
of the interaction but does not strongly infiuence its behavior. The behavior 
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Figure 18. Feynman diagrams that give the pQCD result. Upper panel shows the 
separation of the form factor into the hard scattering part and the initial and final wave 
functions that contain the soft scattering. Lower panel shows one of the thousands (in 
this case) of hard scattering diagrams that make up the hard scattering part. 



is therefore determined by the hard scattering, which can be calculated perturbatively. 
The formalism and method are reviewed in the seminal papers by Brodsky and Farrar 
117|| and Lepage and Brodsky ||118 . 



These calculations of the form factor are all based on the diagrams shown in Fig. 0. 
In the hard scattering, the momentum transfer Q is distributed to the six quarks through 
the five hard gluon exchanges, the last of which carries a momentum of approximately 
Q/6. If the spin factors are included with the quark propagators, the only large 
dependence comes from the l/Q^ of each gluon propagator, giving the counting rule for 
the hard scattering part 

where is the number of constituent quarks (6 for the deuteron) and Uc — 1 the 
number of gluon propagators. This leads immediately to the prediction that the leading 
contribution to the deuteron form factor should go like Q"^^, or that A ~ The 
argument also shows that the perturbative result cannot be expected to set in until 
Q/Q > 0.5 to 1 GeV, somewhere in the region of from 9 to 36 GeV^. (All agree that 
pQCD must give correct predictions at sufficiently high Q^, but how large this must 
be is a topic of considerable controversy ||119| , |12CI|| .) Note that this simple argument does 
not set the scale of the form factor; estimates can be obtained from detailed evaluation 



of more than 300,000 diagrams that contribute to the hard scattering [121] . It turns 



out that this leading twist pQCD estimate is 10^ - 10'^ times smaller than the measured 
deuteron form factor, implying large soft contributions to the form factor, in agreement 
with ||119| , |120|| , and suggesting that pQCD should not be used as an explanation for the 
form factor. The calculation is extremely complicated and a confirmation, or refutation, 
is desirable. 

Perturbative QCD also predicts the spin dependence of the hard scattering, and 
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these predictions provide a more stringent test of the onset of pQCD. These spin 
dependent predictions have imphcations for the individual deuteron form factors, and 
these were first presented in Ref . ||122|| , and further developed in Refs. ||123| , p.24| , |125| . 
Application of these rules to hadronic form factors in general shows that 

• Hadrons with an even number of quarks will be dominated by the longitudinal 
(charge) currents, while these with an odd number of constituents by transverse 
(magnet) currents. Hence, the dominant form factors at large should be the 
nucleon magnetic form factors and the deuteron charge (or quadruple) form factors. 

• The dominant form factors at large are those that conserve helicity. 

When applied to the deuteron, these rules lead to the conclusion that the hehcity 
amplitude Gqq [c.f. Eq. (pO|)] dominates at large Q^; the others are smaller by at least a 
power of Q. In particular, this implies that ~ Gc — 2?7Gq/3 — > {A/Q)Gqq, where 
A is the mass scale above which the nonleading terms can be neglected. Hence 

t2o(Q')^-v^{l + 0(A/g)} . (91) 

Unfortunately, this argument does not allow one to estimate the ratio B/A, since B is 
controlled by a different, independent helicity amplitude. 

In Ref. ||123|| an attempt was made to improve on the constraint (|9lD . These authors 
used the front-form, and evaluated the current in the light-front Breit frame where the 
plus component of the momentum transfer = 0. In this frame all three deuteron form 
factors may be written in terms of matrix elements of the J"*" = J° + component of 
the current 

Gc = {(l - Iv) G+o + IAg+o + i(2r^ - l)G+_} 
Gm = k l2G+o + 2{2r] - 1)^ - 2G'+_| 



G, = nLGt, + 2^-^^GX_Y (92) 



where k = [2p_|_(2?7 + 1)]^^. Perturbative QCD predicts that Gqq will dominate at large 
Q^, and if this happens at a scale A << it follows from (|92D that the form factors go 
in the ratio of Gc '■ Gm '■ Gq = {1 — 2ri/3) : 2 : —1. This leads to a prediction for B/A 
and to a prediction for t2o that differs from ( pT]) at moderate Q^. However, rotational 
invariance is not manifest in the light front, and there are four nonzero components 
of the J"*" current corresponding to deuteron helicity combinations of 00, ++, +0 and 
H — that are related by the angular condition discussed in Sec. |3.6.6| above. Carlson 
has recently shown [ |126| ] that the angular condition places strong constraints on the 
possible subleading behavior of the helicity amplitudes. Perturbative QCD predicts 
that the subleading amplitudes will go like 

GXo - a{A/Q)G+, , G+_ ^ KA/Q^G^o , ^ c{A/Q)'Gto , (93) 

where a, b, and c are dimensionless constants of the order of unity, and A is the scale 
at which pQCD begins working for the deuteron [the Gs in Eq. (^) are identical to the 
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Figure 19. Corrections to the charge form factor, Gc- Each curve shows the difference 
Gcifull) — Gc{nr) expressed as a percentage of the nonrelativistic result Gc{nr). 
The individual cases are discussed in the text. The squares without error bars are the 
extracted Gc data from Ref. and the circles with errors are the Coulomb corrected 
data as discussed in the text. 

Js in Eq. (0)]. Assuming that >> (but making no assumption about the size of 
A) the angular condition in leading order becomes: 

1 + v^a— -f^ G+o = 0. (94) 
\ nid 2mj J 

Solution of this equation therefore requires that A ~ and » m\. Under these 
conditions the relations (p^) again produce only the result (pl|). 

3. 8. Comparison of theory with experiment 

In this section we compare theory with experiment and draw conclusions from this 



comparison. Our major conclusions will be restated and summarized again in Sec. p. 10 
below. 



3.8.1. The charge form factor at very small Figure |T9| shows how coulomb distortion 
of the incoming and outgoing ed plane waves effects the very low data (extracted 
from Ref. [Q) for the charge form factor, Gq- The figure also compares this data with 
theory. 

Corrections for Coulomb distortion change the deuteron radius from an apparant 
2.113 fm (as measured in ed scattering) to 2.130 fm (after the correction) |Q. To 
remove the distortions from the data of Ref. 0, we adjust Gc by ||127|] 

^Gc ^ -0.003 + 0.104 (95) 
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Table 5 



Features of the eight theoretical models reviewed in Sees. 3.6.4, 3.6.6, and 



model 


dynamics 


description 


consistent current 


VOG p9| 


propagator 


bpectator 


yes 


PWM_|| 


propagator 


modified Mandelsweig- Wallace 


no 


FSR |9|~ 


hamiltonian 


instant-form; no v/c expansion 


yes 


ARW_M 


hamiltonian 


instant-form; with v/c expansion 


yes 


CK pof 


hamiltonian 


front-form; dynamical light-front 


no 


LPS liol 


hamiltonian 


front-form; fixed light-front 


no 


AKP JlOQl 


hamiltonian 


point-form 


no 


DB [plSl 


nonrelativistic 


quark-cluster 


yes 



Note that this decreases Gc at very small Q giving a larger deuteron radius, but increases 
Gc where the data have been extracted. The figure shows both the uncorrected and 
the Coulomb corrected data normalized to the nonrelativistic AVIS calculation with the 
MMD nucleon form factor. Note that the difference between the Coulomb corrected and 
uncorrected data is about half of the experimental error at Q — 0.5 GeV. 

The figure also shows the size of relativistic and interaction current corrections that 



arise from the instant-form calculation ARW of Ref. 94 , the front-form calculation 



LPS of Ref. ||104|| , the point form calculation AKP of Ref. ||109|| , and the CIA and 



RIA approximations from Ref. [B9|. These calculations were discussed in Sec. p.6.4| 



and |3.6.6| above. At the scale of the current experimental accuracy (a few percent), 



the relativistic treatments differ noticeably. They also differ from the EFT calculation 
(shown previously in Fig. ^ which drops sharply below the data for Q > 0.2 GeV. It is 
important that these calculations be systematically compared and the different physical 
content of these approaches be isolated and understood. In particular, it would be very 
interesting to know why the covariant CIA and RIA have more positive corrections than 
those obtained from the hamiltonian forms of dynamics. 

3.8.2. Overview of the high predictions The high predictions for the eight models 
reviewed in Sees. |3.6.4|, p.6.6|, and |3.7.1| are shown in Figs. pO|-|2^. The models are 



summarized in Table ^. 

These calculations give very different results. Figure ^ shows the predictions for 
AiQ"^)., with the model dependent p7r7 exchange current intentionally omitted from all 
of the calculations. All of the models except the AKP point-form calculation give a 
reasonable description of A out to ~ 3 GeV^, beyond which they begin to depart 
strongly from each other and the data. Taking into account that the p7r7 exchange 
current could he added to any of these models, and that this contribution tends to increase 
A above ~ 3 GeV^ (for the sign of the p7r7 coupling constant used in the discussion 
in the following paragraph), four models seem to have the right general behavior: the 
VOG, FSR, ARW and the quark model of DB (but there are no results for this model 
beyond = A GeV^). Ironically, none of the models favored by the high Q"^ data does 
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Figure 20. The structure function A for tlie eight models discussed in the text. 
Left panels show the propagator and instant-form results: FSR (solid line), VOG 
in RIA approximation (long dashed line), ARW (medium dashed line), and PWM 
(short dashed line). Right panels show the front-form CK (long dashed line) and LPS 
(short dashed line), the point-form AKP (medium dashed line) and the quark model 
calculation DB (solid line). In every case the calculations have been divided by A 
calculated from the fit (P9) . See Table for references to the data. 



as well at low as the three "unfavored" models shown in the right panels (unless 
the Platchkov H dat a are systematically too low). Another possibility suggested by 
effective field theory is that the asymptotic normalization of the relativistic deuteron 
wave functions is incorrect, and that a small adjustment in NN parameters to insure a 
good value for this constant would correct the problem. 



Figure ^ shows the effect of the new JLab measurements of the nucleon form 
factors on predictions for AIQ"^). These new form factors will decrease predictions for A 
for momentum transfers in a region around = 4 GeV^ (by a factor of 2 for the model 
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Figure 21. Effects of the new JLab measurements and the p-K^ exchange current on 
predictions for A. The sohd hne is the VOG model prediction (in RIA approximation) 
with MMD nucleon form factors (identical to the long dashed curves shown in Fig. |20| ). 
The short dashed line is the VOG prediction (in RIA approximation) with JLab model 
nucleon form factors, Eq. (^). The long dashed curve is the full VOG calculation (CIA, 
/97r7 exchange current with the form factor of Ref. |129|, and MMD form factors). 



shown), increasing the descrepancy between predictions and the data. [However, it may 
improve the prediction for those models (PWM, CK, and LPS) that are currently too 
large in this region.] The figure also shows how the p7r7 exchange current could increase 
predictions at large Q^. The difference between the solid and long dashed lines is due 
largely to the effect of the p7r7 exchange current (but is also do in small part to the fact 
that the CIA result is slightly smaller that the RIA). Unfortunately, the size of the p7r7 
exchange current is very sensitive to the p7r7 form factor, as discussed in Ref. ||128|| , and 
could even be too small to see at these momentum transfers (if the current estimates of 
the p7r7 form factor are too large). This is our reason for insisting that this contribution 
should be viewed as new physics — not readily predictable within a meson model. 

Finally, Fig. ^ shows the predictions for the structure functions A, B, and T20 for 
the eight models discussed. The LPS calculation shows a large descrepency with the 
T20 data, but the most striking feature of these plots is the large model dependence of 
the predictions for B[Q'^). The magnetic structure function provides the most stringent 
test, and the predictions are comparatively free of the p7r7 exchange current (which 
gives only a small contribution to B). Examination of the figure shows that the B 
predictions of the PWM, ARW, AKP, CK models fare the worst. In all, taking the 
predictions for the three structure functions together, the best results are obtained with 
the FSR, VOG, and DB models. 
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Figure 22. The structure functions A, B, and T2o. Tire models, in order of the 
of their mimima in B, are: CK (long dot-dashed line), PWM (dashed double-dotted 
line), AKP (short dot-dashed line), VOG full calculation (as shown in Fig. solid 
line), VOG in RIA (long dashed line), LPS (dotted line), DB (widely spaced dotted 
line), FSR (medium dashed line), and ARW (short dashed line). See Tables]^, |[ and 
^ for references to the data. 
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Figure 23. Upper panels {A and B) and lower right panel (T20) compare data to three 
theoretical models based on VOG: (i) "standard" case referred to in the text (dotted 
line), (ii) model with the tripole F3 (solid line), and (iii) model with the dipole /pir7 
(dashed line). The center two panels show the data and models (ii) and (iii) divided 
by model (i). The lower left panel shows form factors: standard dipole with = 0.71 
(solid line), dipole with = 1.5 (short dashed line), Rome fpT^^ (dot-dashed line) 
[ |l29| , and the tripole with A^ = 5 (long dashed line). 
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3.8.3. Important lessons We conclude our comparison with theory by extending our 
discussion somewhat beyond the hmits of previously published papers. To make the 
point clearly, focus on the VOG calculation using the Spectator equation, and recall 
that current conservation required that the single nucleon current in this approach, 
Eq. (|65D, include a new form factor, F^{Q'^). This form factor must satisfy the constraint 
^3(0) = 1, but is otherwise completely unspecified. In the pubhshed VOG calculation 
p9| and in all of the plots shown so far, this form factor was taken to be the standard 
dipole, Fr,. 

Figure shows how agreement between theory and experiment can be significantly 
improved by choosing a different form for the unknown form factor F3. This figure shows 
three theoretical predictions: (i) the "standard" VOG RIA prediction with F3 = Fd and 
no piT'-f exchange current, (ii) another model with no pTC'j exchange current, but with a 
tripole F3 of the form 

and (iii) a model with F3 = Fd and a dipole p7r7 form factor 

In both form factors is in GeV^. These form factors are shown in the bottom left 
panel of the figure. 



Figure ^ shows that a good agreement between theory and experiment can be 
obtained with the tripole F3 without any p7r7 exchange current, and that to some extent 
the p7r7 exchange current can be substituted for a hard -F3 (although F3 is better than 
/p7r7 at improving all three observables simultaneously). We see that we have, in some 
sense, achived the goals of a theory of elastic ed scattering based on nucleon degrees of 
freedom. With small adjustments of unknown form factors associated with short range 
physics, the NN theory can describe all three form factors quite well. 

It seems likely that any nucleon model with a consistent and complete description 
of the current (c.f. Table ^ can do as well. The reasonable results obtained from the 
FSR and DB models are probably due to the fact that they have consistent, complete 
currents not based on an expansion in powers of (f/c)^ (which must fail at high Q"^). 

3.9. Future prospects 



3.9.1. Future prospects for A From the discussions in Sec. ^]8| above, it is clearly 
of interest to extend measurements of A to higher Q^. An ed coincidence experiment 
is straightforward, but prohibitive timewise with present accelerators. The proposed 
12 GeV JLab upgrade allows one to take advantage of the approximate E"^ scaling of 



crj\/ at constant and high energy ||130|| . A large acceptance spectrometer such as 



MAD would be very helpful. Depending on the details of the upgrade, a one month 
experiment could provide data to of 8 GeV^. 
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It might also be desirable to do a new, high precision experiment at low Q^. The 
goal of this experiment would be to resolve the discrepancy between the data sets 
of Refs. and and to check the low limit of the relativistic calculations. 



These measurements require little time, but do require excellent control of systematic 
uncertainties, at the level of 1 - 2%, if they are to be meaningful. 

3.9.2. Future prospects for B A hasty examination of Fig. ^ might lead one to 
believe that the problems with B are mainly theoretical, and that there is no need for 
new data. We believe this attitude would be inappropriate for two reasons. First, some 
of the calculations shown in the figures are in early stages of development, and will 
improve before any new data are available. The number of calculations is a reflection of 
the challenge that the high data present, and a reflection of interest in, and knowledge 
of relativistic methods that is emerging from the study of these measurements. Second, 
and most important, measurements of B vary by 5 orders of magnitude. The existence 
and position of the first minimum has not yet been firmly established, and the location 
and existence of a possible second minimum is unknown. These minima in B (if they 
exist) result from cancellations of various physical effects and provide a very precise test 
of any theory. 

As indicated above, measurements of B do not require high energy, but do require 
large scattering angles, as close as possible to 180°. At 180°, beam energies of 1 to 2 
GeV cover a range from 1.4 GeV^ to 5 GeV^. 

The SLAG NE4 experiment was a heroic effort, run with ed coincidences. Energy 
resolution was limited by thick targets, 20 - 40 cm long. There was a large but 
manageable background of vr" in the e~ spectrometer, about 3 or 4 to 1, and an 
extremely large proton background from 76? — pn in the deuteron channel - up to two 
events per beam pulse were seen in the worst kinematics. Since NE4 ran at beam currents 
of up to 50 mA instantaneous, with a duty factor of about 0.3 x 10^'^, corresponding 
to an average current of about 15 /iA, the JLab continuous beam structure essentially 
eliminates random coincidence backgrounds such as these. Background coincidence 
reactions included 76? 'yd and 'jd n^d at 180°, with the photon producing an 
electron detected in the electron channel. The spectrum of these photoreactions ends 
near the elastic peak, allowing the background contributions to be fit and determined. 
A JLab experiment would run with both larger spectrometer acceptance and higher 
luminosity to increase rates, but with a shorter target to reduce these backgrounds. 



The proposed configuration ||130|| would use the Hall A septum magnets to detect 
the forward-going deuterons at angles of 3 - 6°, along with special electron channels to 
detect scattered electrons at about 160 - 170°. Based on the SLAG NE4 cross sections, 
a one month experiment can map out B to about 6 GeV^. 



3.9.3. Future prospects for ^20 Extending t2o or other polarization observables to 
higher is quite difficult ||131| , |132|, p.33|| . There are three obvious possibilites. 

Recoil tensor polarimetry requires a well calibrated polarimeter with a large figure 
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of merit, but no such device exists. For example, POLDER, used in the JLab Hall C ^20 
experiment, relied on the ^H[d,pp)n reaction, for which the figure of merit decreases at 
larger energies, leading to a practical upper limit in Q^, that was reached in E94-018. 
HYPOMME ||134| , |135| , p.36|| is promising, but not well enough calibrated. 



The combination of polarized electron beams with recoil vector polarimeters is 
an untested possibility ||137|| . With A, B, and Gm known, pz is calculable and 
can calibrate the polarimeter analyzing power, while Px determines the form factor 
combination Gc + riGq/S. (The ratio of the two polarization components depends on 
this combination of form factors, times kinematic factors and divided by Gm-) Since the 
JLab polarized source can provide 50 - 100 beams, there is no luminosity problem. 
The difficulty with this measurement is that the polarization components are expected 
to be small in the range of interest, ^ 0.01, and Px/Pz ~ 5 - 10. A one month 
measurement for one of ~ 2 - 2.5 GeV^ can determine the polarization components 
well, but if these are small as expected, the form factors can only be extracted with 
factor of two uncertainties. 

An alternative is to use asymmetries from a polarized target. However, the reduced 
currents that can be used with polarized targets require a large acceptance detector 
such as CLAS, to make up for the lack of luminosity. Also, although current polarized 
targets have moderately large deuteron vector polarizations, tensor polarizations are 
small. Furthermore, the asymmetry varies as P2{n-q), with fi the polarization direction, 
so it is desirable to have complete azimuthal coverage, with n in the direction of q 
at one azimuthal angle, rather than being purely transverse. Extensive beam time 
would be needed, either as an external polarized target experiments at JLab, or as an 
internal polarized target experiment at HERMES. We note that a series of moderate 
measurements are planned with the MIT Bates BLAST detector ||138|| , for from 
0.1 to 0.9 GeV2. 

Proposed next generation colliders, such as EPIC/ERHIC, are promising due to 
large planned luminosities; for this experiment the lower proposed cm. collision energies 
are desirable for ensuring exclusivity. The spin direction of the polarized deuteron 
beam must be controllable. In collider kinematics, the scattered electron and deuteron 
energies are close to their respective beam energies and are slow functions of Q^, while 
the scattering angles for fixed vary slowly with the beam energies. Thus, if an 
experiment is possible, it would attempt large azimuthal coverage of coincidence ed 
elastic scattering, with the outgoing particles at angles from a few to about 20° from 
the beam line. 

3.10. Conclusions to Section 3 

Comparison of theory and experiment leads to the following conclusions: 

• Nonrelativistic quantum mechanics (without exchange currents or relativistic 
effects) is ruled out by the ^(Q^) data at high Q^. Reasonable variations in nucleon 
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form factors or uncertainties in the nonrelativistic wave functions cannot remove 
the discrepancies. 

• In some relativistic approaches using A^A^ degrees of freedom only, short range 
physics not calculable within the model (F3 or ff^^, for example) can be adjusted 
to give good agreement with all the data. 

• Different ways of calculating relativistic effects (or meson exchange currents) can 
give results that differ substantially from each other. Even at low Q^, where all 
calculations are constrained, these differences are larger than errors in the data. 
This is not understood, but may be due to the failure of some models to use 
realistic currents. 

• The deuteron form factors provide no evidence for the onset of perturbative QCD, 
but quark cluster models could explain the data. 

Study of the experimental situation leads to the following conclusions: 

• A good database of A, B, and ^20 measurements has been obtained; while 
discrepancies exist they are generally not large enough to affect the theoretical 
interpretation. 

• The minimum of B is very sensitive to details of the models, and improved 
measurements of B for in the region 1.5-4 GeV^ are particularly compelling. 
It is important to accurately map out the zero in the B structure function. 

• Detailed disagreements between theories and different data sets suggests the need 
for precision studies at low Q^. 

4. Deuteron Electrodisintegration 
4-1- Introduction 

As the deuteron has no excited bound states, inelastic scattering experiments have 
largely consisted either of (i) measurements in which the final state mass, W, is very 
close to 2m (referred to as threshold electrodisintegration even when is very large 
because the final state is close to the A^A^ scattering threshold), (ii) measurements near 
the quasifrec peak (defined by the condition that the "spectator" nucleon remain at 
rest), or (iii) deep inelastic scattering in which both and W become very large. 
We will not discuss deep inelastic scattering in this review. For processes at modest 
energies near the quasielastic peak, a rough estimate of the cross section can be made 
using the unrealistic plane wave impulse approximation (PWIA) in which all final state 
interactions are ignored. Denoting the momentum of the outgoing struck nucleon by 
Pi, the cross section in PWIA is proportional to 



where G% is some combination of the squares of the electric and magnetic form factors 
of the nucleon, and (^'(p))^ is an average of the square of the momentum space wave 
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Figure 24. The kinematics of electron scattering when the initial hadronic state is 
broken into two fragments with momenta pi and p2 ■ 

function of the deuteron with internal relative momentum p. Were the PWIA realistic, 
Eq. ( p8D shows that inelastic scattering in quasielastic kinematics would provide a direct 
measure of the (square) of the deuteron wave function. While the PWIA is overly 
simplistic, it does illustrate (correctly) one of the central justifications for quasielastic 
measurements. 

Within the context of a more realistic dynamical theory, one can use response 
function separations and polarization observables to enhance the sensitivity to various 
model dependent nonobservables, such as momentum distributions, meson-exchange 
currents, and medium modifications. One strong recent interest has been to choose 
kinematics in which the unobserved nucleon has a large momentum; the plane wave 
approximation shows that this configuration enhances sensitivity to initial-state short 
range correlations (i.e. the wave function) and possible quark effects. A number of 
these experiments have been carried out at various accelerators, but no experiments 
at JLab have yet reported results. Thus, an experimental review of this topic is 
unwarrented at this time. However, because photodisintegration, electrodisintegration, 
and threshold electrodisintegration are closely related theoretically, we present the 
theoretical background for these processes below. 

4.2. Cross section and polarization observables 

In this subsection the cross section and polarization observables for electrodisintegration 
and photodisintegration to an np final-state are reviewed briefly. (We do not discuss 
pion and meson production.) The electroproduction cross sections will be obtained first, 
and photoproduction will then be treated as a special case. 
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4-2.1. Electrodisintegration The most general decomposition of the d{e,e'p)n 
coincidence cross section was first discussed by Donnelly and Raskin [p,39|| . Here we 
will follow the later work of Ref. ||14(]|] . The cross section can be shown to have the form 
(Eq. (95) of Ref. M ) 

%{r^P + stR? - I \cos 2<PR$1 + sir 
cos (pRix + sin (pRiT + 2/i st'Rt^^ 



dQdE'dE 



6(11) 



+ slt 



+ 2hsLT' sin (pR^l^p, + cos (f)R^l^, 



(99) 



where, in the lab frame, q = {z^, qi}, and the nine response functions are functions of 



Q'^/2mx, and the scattering angle 6i of the final-state proton (measured 



in a coincidence experiment and integrated over in an inclusive measurement). The 
ejectile plane is tilted at an angle with respect to the electron scattering plane, as 



illustrated in Fig. The Mott cross section and other variables are as defined in 
Eq. (Est). The electron kinematic factors are 
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slt' 
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(100) 



with ^ = \qL\Q-Han{9/2). 

One of the virtues of Eq. ( p9D is that the response functions R are covariant, and 
hence (0) it can be used to describe the cross section in either the cm. of the outgoing 
np pair or the laboratory frame, provided we use the appropriate form of (iS: 



dE\ 



dn. 



\LAB = \P1l \ ^ dQi^ , 

where TZ is the lab recoil factor 
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upi — EiqcosOi 



-1 



with W the invariant mass of the outgoing pair. 

In any frame the nine response functions of 

fx) 

the deuteron response tensor R\ y 



(101) 
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are related to the components of 
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All of these quantities are written in the helicity basis, with the helicity of the 
virtual photon, Ai and A2 the helicities of particles 1 and 2 in the final-state, and A^ the 
helicity of the initial deuteron. The matrix element of the helicity-basis current operator 
between helicity states is represented by <^AiA2 J\^{(l) ^dj- In cases where the deuteron 
target might be polarized only in the y direction, and where only the polarization of the 
outgoing particle 1 might be measured, the spin density matrices for particle 1 in the 
final state is given by and that of the deuteron in the initial state by p"', where 

/ 1 -ipf \ 



Nil) 
PXX' 



TV (II) 

Pxy 



V 
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2 







-ipt 
1 



^Py 
1 



{p^{cr.)xx'+P^ia,)xx') (105) 



where Sy is the y component of the spin-one matrix, Py = y3/2ImTii is the vector 
polarization of the deuteron target, Pi is the direction of the polarization of outgoing 
particle 1, measured with respect to the x"y" z" coordinate system shown in Fig. and 
(Ji are the Pauli matrices. Note that only those response functions of type I (denoted 
by the superscript) are nonzero if all of the hadrons are unpolarized; type II response 
functions require (for the cases considered here) measurement of the polarization of the 
outgoing nucleon. Further details and additional cases can be found in Ref . ||140|| . 

The familiar unpolarized inclusive cross section is easily obtained by integrating 
(p9|) in the cm. and summing over electron polarizations. The result is 
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4-2.2. Photodisintegration For real photons the longitudinal components are absent, 
and the cross section simplifies (there is no electron scattering plane and no electron 
kinematics). The most general polarization of the incoming photon, e, is therefore a 
superposition of the circular polarization states e±, which we write as e = a+e+ -|- a_e_ 
with |a-|_p + |a_p = 1. The expansion coefficients a± can therefore be written in terms 
of only three independent parameters 

a+ = -cos/?e-^('^+") 

a_ = sin/Je*^'^-") (108) 
with < /3 < 7r/2. The coincidence cross section for a polarized photon beam is then 
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where p+ = cos 2/? and = sin 2/3 are the fractions of right circular and hnear photon 
polarizations, respectively, and ciSo depends on the frame 

\ / cm. 



dj:o\LAB = dn,^ ( 7^ ) n, (no) 



L 



with TZ defined in Eq. ( |102| ). We will return to the cross section ( |109[ ) in Sec. ^ 

4-2.3. Theoretical Issues It is important to appreciate that these formulae for the 
cross section are exact relativistic results (subject only to the one photon exchange 
approximation) . All of our ignorance is confined to the hadronic matrix elements of the 
current 

^A'\A,(PiP2,g) = (A1A2 |jA,(g)| A,) , (111) 

and the structure functions (|104| ) that are products of these currents. 

In much of the older literature, particularly for studies of the (e, e'p) reaction from 
nuclei with mass number A > 2 [ 141 |, the cross section is written 



d^a 

}Ca,pSip,E,), (112) 



dndE'dJ: 

where /C is a kinematic factor, o"ep is cross section for scattering of an electron from 
an "off-shell" proton, and S{p, Eg) is the proton spectal function (which gives the 
probability of finding a proton with momentum p and separation energy Es in the target 
nucleus). The proton momentum distribution is obtained by integrating the spectral 
function over the separation energy 

n{p) = J dES{p,E) . (113) 

Some early experiements focused on "measuring" the momentum distribution and the 
spectral function. While this picture has a nice physical interpretation [it is motivated 
by the PWIA, Eq. (|98D], and presenting data this way is sometimes useful, particularly 
in the early phases of the program, it is important to realize that the individual structure 
functions that enter the exact cross section ( p9|) are, in general, independent functions 
which are not proportional to each other, and that therefore Eq. (|112D is only an 
approximation to the cross section ||142|| . Attempts to refine the definitions of aep and 



S{p, Eg) can have limited value at best, and at worst can lead to many unproductive 
debates about the precise definition of the spectral function. 

Calculation of the hadronic current matrix elements ( |111|) is complicated by 



requirement that the current be conserved, q^J'^ = 0. For elastic scattering, where the 
initial and final states are identical, invariance under time inversion usually guarantees 
that even simple approximations to the dynamics will satisfy this constraint. But 
building in current conservation for inelastic processes usually requires consistent 
treatment of both final-state interactions and interaction currents. The failure of 
approximate calculations (and the PWIA in particular) to satisfy current conservation 
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is often seen as a serious obstacle. Some ad hoc prescription of the kind introduced by 



De Forest |[L43|| , is needed. 

We propose the simple prescription introduced recently in the study of deep inelastic 
scattering ||144||. Suppose the exact current is composed of two parts = + J2. 



In general, neither of these two parts will satisfy current conservation alone; that is 
Q/xJi 7^ for each i. However, since the exact current satisfies current conservation, 
Qfi{Ji + J2) = 0- We propose replacing each of the individual terms in the current by 

^r-^^. = ^r-^-4^. (114) 

This procedure is covariant, guarantees that each component conserves current (so that 
one can be calculated without knowing the other), and that their sum is unchanged: 



Ji + J2 = Jti + Jt2- Perhaps the best argument can be found in Ref. [T^] where it was 
shown (for a very simple case) that the Born term defined in this way dominates the 
final-state interaction term in the deep inelastic limit, resolving a long standing puzzle. 
Finally, note that Jj^^ = £7^^^ = ^-y^iJri (where are the virtual photon polarization 
vectors satisfying ■ g = 0) so that the response tensor (|103|) is unaffected by the 
redefinition (|114|) ! 

We now turn to a brief discussion of threshold electrodisintegration. 

4.3. Threshold Electrodisintegration 

4.3.1. Overview Threshold deuteron electrodisintegration measures the d{e,e')pn 
reaction in kinematics in which the proton and neutron, rather than remaining bound, 
are unbound with a few MeV of relative kinetic energy in their center of mass system. 
If the final-state energy is low enough, the final state will be dominated by transitions 
to the ^5*0 final state, and will be a pure AS* = 1, AI = 1, Ml transition, similar to the 
— i> A(1232) transition. This transition is a companion to the B structure function; 
both are magnetic transitions and both are filters for exchange currents with only one 
isospin (d — > d is AI = and d —>■ ^Sq is AI = 1). To see the similarity, compare the 



top right panel of Fig. ^ with the threshold measurements shown in Fig. Both have 
a similar shape, and in both cases the uncertainties in the theoretical predicitions are 
large. 

The similarity of these two processes (elastic and threshold inelastic) also holds 
for the theory. These two processes can be used to separately determine the precise 
details of the 1 = and 7 = 1 exchange currents. Once the exchange currents are fixed, 
they can be used to predict the results of d{e, e'p)n over a wide kinematic region. Any 
theoretical approach that works for the form factors should also work equally well for 
threshold electrodisintegration, yet very few of the groups who have calculated form 
factors have also calculated the threshold process. This may be due in part to the 
fact that the 7 = 1 interaction currents are larger than the 7 = interaction currents, 
making the threshold electrodisintegration calculation more difficult than the elastic 
calculation. A more definitive test of the various relativistic approaches discussed in 
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the previous sections will be possible once the elastic calculations are extended to the 
threshold inelastic process. 

Previous threshold electrodisintegration experiments have reported an average cross 
section that can be obtained theoretically by integrating the relative pn energy in the 
final state, Epn, from to 3, or (in some cases) to 10 MeV. The unbound ^5*0 final 
state dominates at threshold (because the ^Si-^Di scattering state is orthogonal to the 
deuteron state at threshold), but above threshold there are contributions from the ^5*1- 
^Di scattering state, and eventually from the A^A^ L = 1 scattering states as well. To 
emphasize this magnetic transition, data have been taken at large electron scattering 
angles, mostly 155° or 180°. 

Threshold electrodisintegration provides strong evidence (perhaps the best we have) 
for the existence of isovector exchange currents |145|. The impulse approximation 



calculation of the transition to the ^Sq final state has a zero arising from the negative 
interference between the ^Si Sq and ^Di — 5*0 pieces of the transition that lead to 
a minimum at near 0.5 GeV^. This minimum is not seen in the data, and theoretical 
calculations of the 1 = 1 exchange current contribution fill in the minimum and explain 
the data. 



4-3.2. Current status of theory and data Experiments have measured d{e,e')pn 
threshold electrodisintegration to about 1.6 GeV^ with better than 1 MeV (cr) 
resolution, integrating up to Epn = 3 MeV |146| , p.47| , |148| , |149| , p.5U| |, and to 



nearly 3 Gev^ with 12 - 20 MeV (FWHM) resolution, and integrating up to Epn = 10 



MeV [|m|, ^ 



Figure ^ shows the smooth rapid fall off of threshold electrodisintegration cross 
sections for Q"^ up to about 1.2 GeV^, and the change in slope for the higher SLAG 
and Bates data. There is good agreement between the various measurements, including 
those not shown, considering the change in scattering angle (155° for Saclay, 160° for 
Bates, 180° for SLAG) and integration region for Epn (0-3 MeV for Saclay and Bates 
data shown, 0-10 MeV for SLAG). The figure also shows results from seven theoretical 
calculations. 

This model 



The oldest calculation shown in the figure is from Mathiot |153 



includes one and two body currents based on vr and p exchange, and also contributions 
to the two-body current operator from the exitation of the A. Mathiot confirms that 

II One possible source of confusion in any close examination of threshold electrodisintegration is the 



use of two different conventions for the cross section. Early articles |147, 148|] and the most recent Bates 
article [ |15C| report the cross section as (Pa/dfldiv = d^a/dfldE' where u {E') is the energy transfer 
(scattered electron energy). The SLAG articles |151, |152| and first Bates article |149| use instead 
d^a/dfldEnp, where Enp is the total np kinetic energy in their cm. frame, Enp = W—md — 2.225(MeV). 



(A typographical error at one point in |152| misidentifies Enp as the energy of a nucleon, rather than the 
two nucleons.) These cross sections are related by the Jacobian \dEnp/dE'\ = {rrid + 2E sin^ {6 /2)) /W 
which numerically ranges from about 1.3 to 2.3 for the data that we present. The articles showing the 
"dE'np" cross sections appear to plot the Saclay cross sections as published, as "di?"' cross sections, 
rather than converted. 
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Figure 25. The cross section for threshold electrodisintegration of the deuteron. Data 
shown are from Saclay B [p7| , Saclay A [psl , Bates |l49i [l5C|] , and SLAG [Tsil, 



Theory calculations are from Mathiot |153| , LSA 
YYW and LC |58|. 



SchiaviUa |55|, STG |15§ 



the one body current (the impulse approximation) produces a sharp minimum in the 
cross section at about = 0.5 GeV^ (in agreement with the recent Schiavilla-1 body 
curve shown in the figure and in complete disagreement with the data) and that the vr 
exchange current fills this in, shifting the minimum to about 1 GeV^. The p exchange 
is also important, shifting the mimimum to ~ 1.4 GeV^. Contributions from the 
electroexication of a A are smaller, at least below — 1 GeV^. Figure ^ shows 
that this calculation breaks down above ~ 1 GeV^, probably because it does 
not include many of the contributions included in modern calculations. More recent 
calculations, also based on hadronic degrees of freedom, are from Leidemann, Schmidtt, 
and Arenhovel [LSA] ||154|| and Schiavilla. (The calculations by Schiavilla shown on 
the figure are based on the work of Ref. but use the more recent Argonne AV18 
potential. These curves were also published in the review ||155|| .) In both of these 
cases the details are largely unpublished and the inelastic calculations are not as up- 
to-date as the corresponding elastic calculations recently done by the same groups, so 
it is premature to draw definite conclusions. The work of Smejkal, Truhlik, and Goller 
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[STG] ||156|| obtains exchange current contributions from the vrpoi system using a chiral 
lagrangian, and does a good job describing the data out to 2 GeV^. Predictions from 
two quark cluster models, the early model of Yamauchi, Yamamoto, and Wakamatsu 
[YYW] ||157|] , and the more recent model of Lu and Cheng [LC] ||158|| , are also shown 
in Fig. These calculations both tend to have too much structure in the region of 
the shoulder, but do show that quark cluster models have the ability to describe the 
exchange currents needed to account for the data. 

The recent improvements in relativistic theory discussed in Sec. |3.6| will lead to a 
new generation of calculations that will rely on threshold electrodisintegration to provide 
details about the nature of the 1=1 exchange currents. The most precise constraint 
on these currents comes from the d ^ ^Sq transition, and this part of the transition is 
partly obscured by the poor energy resolution of the existing high measurements. A 
new and improved experiment at JLab with higher resolution would allow the threshold 
d — > process to be better extracted, with a better resulting determination of the 
isovector exchange currents. It is also important to determine whether or not there is 
a minimum near 1.2 GeV^. This present indication of a minimum might be an artefact 
of the end of the Saclay data ||148|| vs. the start of the SLAG data ||151| , |152|| , along with 



systematic uncertainties of these and the Bates measurements ||149| , p.5CI| . 

A high statistics high resolution measurement at large is feasible. Measurements 
have been proposed ||159|| with 1.5 MeV resolution at a scattering angle of 160°. The 
experiment would use Hall A with the HRS spectrometer vacuum coupled to the 
scattering chamber, and a special cryotarget and collimators to enhance resolution and 
reduce backgrounds. Measurements were proposed for 6 points from 1 to 3.7 GeV^. 

We note that SLAG NE4 simultaneously measured the threshold electrodisintegra- 
tion along with the elastic structure function B. Because the large threshold inelastic 
cross section is typically an order of magnitude larger, the d{e, e') measurements see es- 
sentially only the inelastic processes, and d{e, e'd) is needed to determine the elastic 
scattering. If it is possible to maintain large solid angle for the elastic scattering, and 
good resolution for the threshold electrodisintegration, both data sets can be obtained 
simultaneously. 



5. Deuteron Photodisintegration 
5.1. Introduction 

Deuteron photodisintegration was first investigated in the early 1930s, in order to 
understand the structure of the neutron. After the discovery of the neutron by James 
Ghadwick, attention turned to its mass and structure. Was the neutron a fundamental 
particle, like the proton and electron, or was it a bound state of the electron and 
proton, different from the hydrogen atom? If it was a bound state of the proton and 
electron, how were the electrons confined into the small nuclear volume? Gonfiicting 
experimental evidence on the neutron mass prevented resolution of the issue until 1934, 
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when Chadwick and Maurice Goldhaber used deuteron photodisintegration |p.6CI|| to 
determine that the neutron mass was shghtly heavier than that of the hydrogen atom. 
Thus, the neutron, being heavier than the proton plus electron, was a fundamental 
particle, and there was no longer any basis for thinking electrons could be present in 
nuclei [p.61| . 

Subsequently, deuteron photodisintegration cross sections have served as a standard 
test case for nuclear theory. The effects, for example, of the A resonance in cross 
sections for beam energies near 300 MeV are pronounced, but until recently large 
discrepancies between different experimental data sets made precise tests of theories 
difficult ||162|| . The 1960s-1970s saw the start of polarization measurements. The earliest 
data were intermediate energy measurements of the induced proton polarization [|163 



and low energy measurements of the induced neutron polarization ||164|| . Large induced 



polarizations were observed [|165| , |166|| soon afterward, particularly for energies above 



the A resonance and for center of mass angles near 90°. The combination of more 



extensive confirming measurements ||167|| for Ej about 350 - 700 MeV, which could 
not be reproduced theoretically, and interest in dibaryons led to much excitement about 
deuteron photodisintegration in the late 1970s and early 1980s. There were many serious 
theoretical efforts, numerological studies involving inclusion of dibaryon resonances, and 
extensive experimental studies of cross sections and polarization observables. 

5.2. Relation between elastic scattering and photodisintegration 

Since both the recent deuteron form factor measurements and the recent high energy 
deuteron photodisintegration measurements have been made with 4 GeV electron beams, 
it is sometimes assumed that the same theory should work for both. In this review 
we emphasize that this need not be the case. The kinematics of elastic electron- 
deuteron scattering and deuteron photodisintegration are very different, and the physics 
being explored by these two measurements is also very different. The implications of 
this remarkable feature of electronuclear physics, often not fully appreciated, will be 
discussed briefly in this section. 



The kinematics of elastic scattering and photodisintegration are compared in Fig.p6 
which shows W function of the photon (real or virtual) energy 

W^-ml = 2mdU (l - — ) , (115) 
V mdJ 

where x = Q'^/2mv and v = E^ for real photons. The mass of the final excited state 
increases rapidly as x decreases below its maximum allowed value of a; = m^/m ^ 2. 
For any energy u or any Q^, elastic ed scattering leaves the pn system bound, with no 
internal excitation energy added to the two nucleons. For quasifree scattering (x = 1) 
the mass of the final pn system grows with i/, and as x decreases below 1 the mass 
grows more rapidly with v. As x 0, we approach the real photon limit. Real photons 
produce the maximum value of W of any given beam energy. With each 1 GeV of 
beam energy W — rud increases by approximately 500 MeV, driving the final state deep 
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Figure 26. The variation of with the photon energy v for various values of 



X, 



as given by Eq. (115). The shaded regions show the approximate thresholds for the 
production of bands of nucleon resonances, as discussed in the text and shown in Table 
^ The numbers in the small circles are the number of distinct channels in each band. 

Table 6. The 24 well established nucleon resonances listed in the Particle Physics 



Booklet 1 168 1 fall into the 8 bands listed below. Masses of neighboring resonances in 
each band are less than 150 MeV apart. All of these resonances can contribute to 
deuteron photodisintegration for W < 4.5 GeV. All but A^4 and A4 can contribute in 
all combinations, giving 13 + (13 x 12)/2 — 91 channels with two 1 — 1/2 particles, 
45 channels with two / — 3/2 particles, and 117 channels with one 1=1/2 and one 
7 = 3/2 particle. The number of additional channels contributed by A'4 and A4 is 
shown on the table and totals 33 channels. The total number of channels is 286. 





I = 1/2 




7 = 3/2 




N (939) 


Ai 


7^33(1232) 




Fii(1440), 7^13(1520), 5ii(1535), 


A2 


P33(1600), ^31(1620), 7^33(1700) 




5ii(1650), 7^15(1675), Fi5(1680), 








7^13(1700), Pii(1710), Pi3(1720) 






m 


Gi7(2190), 77i9(2220), Gi9(2250) 


A3 


F35(1905),P3i(1910), P33(1920), 








7^35(1930), F37(1950) 




7i4i(2600) 14 channels 


A4 


73,11(2420) 19 channels 



into the resonance region. The well established nucleon resonances, all of which can be 
excited by 4 GeV photons, are listed in Table ^ and the bands of thresholds at which 
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Table 7. The thresholds for the production of pairs of baryon resonances also fall 
into 8 bands. Neighboring thresholds within each band are less than 150 MeV apart. 



These bands 


are shown in Fig 


HI- 




Band 


1.1J.CLOD i- dJ-itiC 


TTIPITI nPTQ 


■nninnPT' of 










Bi = NN 


1878 


NiNi 


1 




2171 


iViAi 


1 




Z404 — ZO (i) 


A A AT AT AT A 
AiAi, i\'2iV2, JV1A2 




B4 


2858 - 2872 


7V2A1, N1A3, A1A2 


17 


B5 


3155 - 3280 


A1A3, N2N2, N2A2, A2A2 


86 


Be 


3452 - 3652 


7V3A1, 7V1A4, iVliV4, ^2A3, 








A2A3, A1A4 


66 


Br 


3832 - 3860 


7V4A1, A3 A3, N2N3, N3A2, 


52 


Bs 


4046 - 4440 


7V2A4, A2A4, iV3A3, N2Ni, 








A2N4, A3A4, N3N3 


50 



these resonances are excited, either singly or in pairs, are hsted in Table |^ and shown in 
Fig. |26|. [The Fermi momentum of the struck nucleon, and the widths of the resonances, 
will average these thresholds over a wider kinematic region than shown.] By = 1.2 
GeV, the final-state mass is already reaching W = 2 GeV, the nominal onset of deep 
inelastic scattering (DIS) from a single nucleon if we assume one of the nucleons in the 
deuteron remains at rest. At E^ = A GeV, the final-state mass is approximately 4.5 GeV, 
and at least 286 thresholds for the production of pairs of baryon resonances have been 
crossed (and there are probably more from unseen or weakly established resonances). A 
photon energy of 4 GeV corresponds to np scattering with a laboratory kinetic energy 
of about 8 GeV! (See also the kinematic argument given by Holt ||169|| .) 



It is clearly very difficult (if not impossible) to construct a theory of high 
energy photoproduction in which all of the 24 established baryon resonances and their 
corresponding 286 production thresholds are treated microscopically. By contrast, 
elastic electron deuteron scattering requires a microscopic treatment of only one channel. 
All of the 286 channels also contribute to elastic scattering, of course, but in this 
case they are not explicitly excited, and can probably be well described by slowly 
varying short-range terms included in a meson exchange (or potential) model. In 
photodisintegration, each of these channels is excited explicitly. As we shall see below, it 
is sufficiently difficult to construct an adequate theory in the region of the A resonance, 
so it is difficult to imagine this program being extended to a realistic treatment of 
many resonances, including intermediate states in which multiple mesons are present. 
It would instead appear that an alternate framework that averages over the effects of 
many hadronic states is needed. The alternatives are to use a Glauber-like approach, 
or to borrow from our knowledge of DIS and build models that rely on the underlying 
quark degrees of freedom. We will return to these issues in our review of the theory in 
Sec. |5^ below. 
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5.3. Cross section and polarization observables 

The photodisintegration cross section was previously given in Sec. [4.2. 2| . Including the 
polarization observables defined by the density matrices of Eq. ( |105| ), and using the 
notation of Ref . ||140[ , the structure functions are 



d(I) 



R? «, IniTn, U) = Rt + Rriv) + pt RrilmT, 



y 

?(i) I^N 



III 



?(II) _ 6(11) /„JV 



Rrp, = Rj'l {p^l^P^I 



R. 



(11) 

TT' 



5(11) ( N 
J^TT' \Px' ' Pz' ) 



RxT + Py Rxriy) 

= p^, Rt{x') + pf, Rt{z') 
= P^,RtAx')+p^,Rtt{z') 



(116) 



where the coordinate system, shown in Fig. (a simplified and relabeled version of 
that given in Fig. |2^) is constructed from the incident photon direction k = z and 
the outgoing proton direction p = z'. Substituting the expansions (|116|) into the cross 
section formula ( p.09| ), gives, in a notation suggested by Ref. |p.70|| . 



da 

dn 




N 

PyPy 



PyT - p^{J: + PyTi) cos 20 
+P+{Cx'Px' + C^:p^>) +PjiO.^>p,^> + 0^'P^ 



sin 20 



;ii7) 



where {da/dfl)Q is the differential cross section for unpolarized photons, and explicit 
expressions for the asymmetry parameters are given in Table |[ Note that the 
observables Cx and Py are the real and imaginary parts of the same combination of 
amplitudes, so that an experimental measurement of both of these observables will fully 
determine this linear combination of amplitudes.0]] This is not true for any other pair of 
observables shown in the table. In the cm. system 

/d(T\ 

[dnj. 



8uoW 



R 



a Pi 



m 



:il8) 



where 



(119) 



The Fs are defined in Table |^ 

The single polarization observables are the induced proton polarization py, the 
linearly polarized photon asymmetry S, and the vector polarized (along y) target 
asymmetry T. The quantities Py and E are defined by 



Py 



/ da/dVL^ 
\da /dVt^ 



da/dQ^ 



da/dVt, 



1 I da / dVt — da / dVL\\\ 
+ da/dVt\\j 



p^ \da/dVL 



(120) 



where ± refers to = ±1, and in the expression for S the photon is polarized either 
in plane along the x direction (0=0, denoted ||) or out of plane along the y direction 
(0 = 7r/2, denoted ±). Note that the opposite phase convention for {a\\ - a<^) is often 
used, and is used in the figures below. The double polarization transfer observables 



^ Some references, such as [170 , use Py for the induced polarization instead of Py. 
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Figure 27. The cordinate system for deuteron photodisintegration. The y = y' axis 
is given hy y = k x p/\k x p\, and the x, x' axes are chosen to make a right handed 
coordinate system. 

Table 8. Formulae for polarization observables. Each structure function in the 
second column is to be divided by Rt, and in the third column by f{0). [The Fi± are 
defined in Table 

Observable Structure Helicity amplitude combination 
function 

Py Rriy) 
T i?T(ImTii) 

Ti Rrriy) 

C7,, Rt{x') 

C,, Rt{z') 

O^i Rtt{x') 

O.' Rtt{z') 



are Ti and Ox',z' (for linearly polarized photons), and C^'^z' (for circularly polarized 
photons), with the subscripts on O and C giving the polarization direction in the final 
state. 

5.4- Theory 

5.4-1 ■ Models using meson-haryon degrees of freedom There have been many attempts 
to understand low energy deuteron photodisintegration using a conventional meson- 
baryon framework. Since the first band of nucleon resonances is not excited until about 



2ImE 



i=l 



(i+3)+ 



2ImELi E]=o [F(i+3j)+ F*i+3j+i)+ + F{i+3j)- -^(*+3j+i)- 



2ReE?=i(-) 
2ImELi(-) 

ELA\F^;\'-\ 

2ImELi(-)^+^ 
2ImELi(-)^+^ 



-Fi+ F*^_iy 



+ F, 



{3+j)+ 



-F- , F* -4- F- F* 



FU Fii+3y 



+ Fi~ ^(*+3)+ 



Z? 77* I Z? 77* 



Fi+F, 



(4-i)- 



+ F(3+i)+ Flj_.i)_ 
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Table 9. The relations between the hehcity amphtudes used in Refs. [140| and |17C] 



Amphtudc 



Ref. 1 140 1 



Ref. |17C| 



(±i±i 


1-^ 


e+ 


|1) 




Fi±/2to 


(±i±i 


1-^ 


e+ 


10) 


-^3,4 


F2±/2to 


(±i±i 




e+ 


1-1) 


-^5,6 


F3±/2to 






e+ 


|1) 


-^7,8 


F4±/2to 




1-^ 


e+ 


10) 


-^9,10 


F5±/2m 




1-^ 


e+ 


1-1) 


-^11,12 


F6±/2to 



400 MeV photon energy (recall Fig. ^) it makes sense to describe the process below 
400 MeV using a model of coupled A^A^, A^A and NNn channels. Laget ||171|| showed 
that the prominant shoulder in the total cross section at = 300 MeV can be largely 
explained by the mechanism in which a A is photoproduced at a nucleon followed by the 
reabsorption of its decaying pion by the other nucleon. He also examined many other 
mechanisms, including rescattering up to second order, but did not do a full calculation 
of the final-state interaction. Later Leideman and Arenhovel 111 7211 treated the A^A^, 



A^A and AA as coupled channels and included final-state effects to all orders. Tanabe 



and Ohta ||1 73|| followed with a more complete treatment of the final state which is 
consistent with three-body unitarity. In a number of conference talks, Lee [|174|] reported 
on coupled channel calculations using A^, A and the Pii(1440) (Roper) resonances, which 
he extended to ii^-y = 2 GeV. His work suggests that final-state interactions significantly 
enhance the cross section for photon energies above 1 GeV. 



The recent calculations by Schwamb and Arenhovel and collaborators ||175|| include 
the A^A^, A^A and nd channels, and also contributions from meson retardation, meson 
exchange currents, and the meson dressing of the nucleon lines required by unitarity. 
All parameters are fixed from nucleon-nucleon scattering and photoreactions such as A 
excitation from the nucleon, so no new parameters are introduced into the calculations of 
the deuteron photodisintegration process itself. They obtain a reasonable description of 
A^A^ scattering up to lab energies of 800 MeV, particularly for the important partial 
wave, and emphasize that the consistent inclusion of retardation effects improves their 
results for photodisintegration. In another work it was shown that inclusion of the 
5*11 and Di^ resonances | |176| | seems to have only a small effect below about 400 MeV. 
These resonances enhance the total cross section by only 14% at 680 MeV, although 
effects on double polarization observables can be more significant. Hence there is some 
justification to limiting the theory of low energy photodisintegration to the channels 
considered in Ref. 



175 



By comparison to the detailed and careful treatments developed for lower energies, 
calculations specifically designed to describe higher energy photodisintegration are less 
complete. An unpublished Bonn calculation [|1771| includes pole diagrams generated 
from TT, p, rj, and u exchange, plus the 17 well-established nucleon and A resonances 
with mass less than 2 GeV and J < 5/2 (listed in Table |^). The calculation uses 



Electromagnetic structure of the deuteron 



71 



resonance parameters taken from the Particle Data Group. Nagornyi and collaborators 



178, 179 1 have introduced a covariant model based on the sum of pole diagrams. The 



latest version ||179|| gives the photodisintegration amplitude as a sum of contributions 



from only 4 Feynman diagrams: three pole diagrams coming from the coupling of the 
photon to the three external legs of covariant dnp vertex plus a contact interaction 
designed to maintain gauge invariance. This model divides the dnp vertex into "soft" 
and "hard" parts, with the hard part designed to reproduce the pQCD counting 
rules and its strength determined by a fit to the data at 1 GeV. The model has no 
final-state interactions or explicit nucleon resonance contributions. There is also a 
relativistic calculation of photodisintgration in Born approximation using the Bethe- 
Salpeter formalism [|180| , |181|| . It is found that the cross section is a factor of 2 to 10 times 



too small, and that relativistic effects are large. All of these models are rather crude, 
and taken together it is not clear what one should conclude from them. The calculations 
each appear to emphasize some aspects of what a comprehensive meson-baryon theory 
of photodisintegration should entail. Perhaps we can say that conventional calculations 
that neglect final-state interactions seriously underestimate the cross section, but may 
be corrected in an approximate manner by introducing diagrams with poles in the s 
channel. S-channel pole diagrams can be regarded as a crude approximation to the 
missing final-state interactions. More generally, using pole diagrams with form factors 
that have the correct behavior at high momentum transfer may also insure that meson- 
baryon theories of deuteron photodisintegration will also have the correct high energy 



behavior 1182 



5.4-2. Models based on quark degrees of freedom Perturbative QCD, discussed briefly in 
Sec. p.7.2| , provides explicit, testable predictions for the cross section and polarization 
observables. For the case of elastic scattering, the high energy (virtual) photon had 
to share its mommentum equally with all of the constitutents, leading to the typical 
diagram shown in Fig. 0. Photo (or electro) disintegration differs in that the momentum 
will not be distributed to all six quarks unless momentum is transferred to each nucleon, 
and this requires non-forward scattering, i.e. the angle ^cm between the three-momentum 
of the outgoing proton and the photon (in the cm. system) must not be or vr. More 
precisely, if is the energy of the photon in the lab system, the square of the momentum 
transferred to each nucleon (where i = 1,2) can be written 

t= {q- pif = m^ - TUdE^ \ 1 - y^l - cos 6'cm j -s sin^ \9cm 

t2 = u = {q- = m^ - nidE^ ( 1 + ^1 - cos 6'cm j ^ s cos^ \9cm , (121) 

where s, the square of the cm. energy, is 

s={d + qf =ml + 2mdE^, (122) 
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Figure 28. A "typical" pQCD diagram describing large angle photodisintegration. 
There are 5 hard gluons. The RNA model places the hard gluons in the rectangular 
boxes within nucleon form factors. 



and the limits in ( |121|) are the resuh for ^ m^. If ^cm = 90° the momentum 
transferred to each nucleon is balanced and the pQCD result is reached most rapidly. 
A "typical" pQCD diagram leading to large angle scattering is illustrated in Fig. In 
general we cannot expect the cross section to follow pQCD unless the minimum of the 
two momentum transfers ti and ^2 is larger than some value tmin at which pQCD holds. 

If the above conditions are satisfied, then the cross section and polarization 
observables should satisfy results predicted by pQCD: the constituent counting rules 
(CCR) and hadron helicity conservation (HHC). 

The CCR | |117| , |183[ predicts the energy dependence of scattering cross sections at 
fixed center of mass angle 

| = ^/(M. (123) 

where n is the total number of pointlike particles in the initial and final states of the 
reaction. For deuteron photodisintegration, n = 13 (there is one photon and 6+6=12 
quarks), and because of Eq. ( |121[ ) t may be replaced by s, as is commonly done. For 
elastic ed scattering n = 14 because there are electrons in both the initial and final 
state, and t = Q^, the momentum transfered by the electron. In this case dt ~ Q^dQ, 
and recalling Eq. ( PB| ) we have 

/ . \ 12 



(124) 



recovering the pQCD prediction A ~ Q^^^. 

The second rule from pQCD is that the total helicity of the incoming and outgoing 
hadrons should be conserved ||118| 

(125) 



where Aj (A/) are the helicities of the initial (final) hadrons. In Sec. |3.7.2| we discussed 
the implication of this rule for the deuteron form factors. HHC makes predictions for 
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Table 10. Implications of HHC for some polarization observables in 70? — > pn. The 
90° limits require additional assumptions. 



Observable HHC limit Approach 
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many spin observables, particularly vector polarizations ||184|| . Here it predicts that the 
amplitudes of Table ^ will have the following behavior at large t 

F3_, leading 
F2-1-, F4-1-, Fq± suppressedby t^^ (126) 

F3+ suppressedby 

where we have used the fact that each hehcity flip is suppressed by a power of t. The 
implication of HHC for the observables given in Table ^ is summarized in Table 
The limits at O^ra = 90° on some of the observables require assumptions about relations 
between the helicity conserving variables ||178| , |185|| . 

In an attempt to include some of the soft physics omitted from pQCD, and to 
extend the region of applicability of pQCD down to lower momentum transfers, Brodsky 
and Hiller introduced the idea of reduced nuclear amplitudes (RNA) ||186|| . In this 
model the gluon exchanges that contribute to identifiable subprocesses (such as nucleon 
form factors) are collected together and their contributions replaced by experimentally 
determined nucleon form factors. It is hoped that the resulting expressions will correctly 
include much of the missing soft physics, and will therefore be valid to lower momentum 
transfers than the original pQCD expressions from which they were obtained. When 
applied to deuteron photodisintegration, the cross section is written 

f = 2A7r^i7-miy ^ ^ (7^^ ^^'^'"^^'^^"^^^'^^'^"'^ ^^^^^ 
where the phase space factor of l/(s — m'^Y comes for a careful reduction of the phase 
space factors in Eqs. ( |103| ), ( |110| ), and ( |121| ) f{6cm) is the reduced nuclear amplitude, 
Fp and F„ are the proton and neutron form factors with tp and t„ the average momentum 
transferred to the proton and neutron 

ip = {pi-d/2f i^ = {p2~d/2f, (128) 
+ In some experimental fits the phase space factor used is (s — m^)^/'^\/ s — Am? instead of (s — m^)^. 
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sin 9c 



(129) 



and the square of the transverse momentum is 
' s 
A 

The power of the "extra" factor of is fixed once the phase space and nucleon 
form factors have been taken into account. Note that this model does not attempt 
to normahze /. 

The RNA form given in Eq. ( |127| ) is somewhat arbitrary, particularly in the 
specification of the form of p"^ Q. Radyushkin ||187|| has argued that the elementary 
process not accounted for by the nucleon form factors (i.e. the absorption of a hard 
photon followed by exchange of a hard gluon with another quark) should include 
nonperturbative contributions. The effect is to replace the l/p^ factor in Eq. ( |127| ) 
by a smooth function /^(s, t), which is assumed to vary slowly in energy and angle. 
In the fits described below, p will be taken to be a constant adjusted to fit the data, 
implying that da/dt ~ s~^° instead of s~^^. 



Alternatively, if the quark exchange mechanism shown in Fig. ^ is to be taken 
seriously, a more detailed calculation is possible. This is the motivation for the work of 
Frankfurt, Miller, Sargsian, and Strikman ||188| , |189|| where the quark exchange diagram 
(which the authors refer to as a quark rescattering diagram) is calculated in front- 
form dynamics using model wave functions for the nucleons and the deuteron. The 
matrix element is written as a convolution of an elementary quark exchange interaction 
with the initial and final nucleon wave functions. The final nucleons are free and the 
distribution of the initial nucleons is given by the deuteron wave function. Since the 
photon momentum is shared by the proton and neutron, there is little sensitivity to the 
high momentum part of the deuteron wave function. The elementary interaction is a 
quark exchange between the two nucleons, with the photon absorbed by one quark which 
then gives up its momentum through a hard gluon exchange with another quark. The 
authors show that this can be replaced approximately by the wide angle np scattering 
cross section (also dominated by quark exchange). The final formula obtained from 
these arguments is 



da 
It 



■yd 



9(s 



■C [ip/s 



da{s, 



dt 



np 



{2tt 



(130) 



where C{—x) = x/{l — x) was used in the fits to the data discussed below. The authors 
propose using experimental data for the np cross section, but since data does not exist 
for the actual kinematic conditions needed, it must be extrapolated, and predictions for 
photodisintegration are given as a band corresponding to the uncertainties introduced 
by the extrapolations. The authors believe that their predictions should be valid for 
> 2.5 GeV, and nucleon momentum transfers —t = —ti and —u = —t2 > 2 GeV^. 

The high energy approaches described above all focus on the region where both 
ip and ip are large (where perturbative arguments can serve as the foundation for the 



* In Ref. [186| it is proposed that p'^ = tu/s. This choice of pt leads to an increased energy dependence 
and worse agreement with the data. 
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treatment). Alternatively, we may ask what to expect when one of these momentum 
transfers is small (but s is still large). The authors in Refs. ||190| , |191|| develop a model 



(131) 



(which they refer to as the "quark-gluon string model" ) based on a Reggie generalization 
of the nucleon exchange Born term. Here the exchanged nucleon is replaced by a nucleon 
Reggie trajectory that represents the sum of a tower of exchanged nucleon resonances 
(or the exchange of three quarks dressed by an arbitrary number of gluons). The energy 
dependence of the predicted cross section is 
da /s\2Qjv(t)-2 

dt V So / 

where aN(t) is the nucleon Reggie trajectory, with aAr(t) = — 0.5 + 0. 9t + 0.25t^/2, where 
t is in units of GeV^. Recent work [|191|| emphasizes the importance of the nonlinear 
term in the Regge trajectory. The model is intended to work for > 1 GeV, with —t 
less than about 1 GeV^. 

We now turn to a review of the experimental data and to a comparison between 
theory and experiment. 

5.5. Experimental Status 

The world data set for cross sections (except for the most recent experiments) has been 



presented several times [162], and we will not review the normalization problems of older 
cross section data sets. Table |ll] presents an extensive list of the published polarization 
data. About 70 pubhcations, starting in 1960, present about 1200 data points for 



photodisintegration and the time reversed radiative capture reaction. Table [11] generally 
lists photon lab energy and proton cm. angle (neutron cm. angle for the Py data). All 
of the radiative capture experiments have measured Ay, and for these experiments we 
usually give the neutron beam kinetic energy and outging photon cm. angle; for the 
lUCF experiment we give their reported proton cm. angle. Matching cm. energies leads 
to T„ = E^md/rnp+{m'^ — {mn+mp)'^)/2mp or T„ ^ 2E^ at high energies. Comparison of 
these data indicates serious problems with backgrounds and/or estimates of systematic 
uncertainties in a number of cases, as will become clear in figures in the sections below. 

We review some of the lower energy data in the next subsection. High energy 
experiments, with photon energy above ~1 GeV, are covered in the following subsection. 

5.5.1. Low energy photodisintegration In this section we review selected experiments 
with beam energies from about pion production threshold to several hundred MeV. 
Tagged photon facilities, with their improved knowledge of incident beam flux, have 
allowed significantly improved cross section measurements in this region since the 1980s. 
Of particular note are extensive recent data sets from LEGS and Mainz. Figures p9|and 
pO| show angular distributions for the cross section, and S, T, and Py at photon energies 
near 300 and 450 MeV, respectively. 

The LEGS tagged, backscattered, and linearly-polarized photon beam was used 
to determined cross sections and S ||248| , |249| , |250|| . Five independent measurements 
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Table 11. Measurements of deuteron photodisintegration polarization observables. 



Laboratory 


Observable 


tj (^iVleVj 


"cm V ) 


^ of Points 


MIT 


Py 


250 


49 


1 


Livermore 


Py 


2.75 


50 - 136 


5 


Zurich 


Py 


2.75 


44, 94 


2 


Illinois 


Py 


12 - 23 


148 


4 


RPI 


Py 


12 - 30 


90 


3 


Purdue 


Py 


294 


72 


1 


Stanford 


S 


80 - 140 


45, 90, 135 


41 


Livermore 


Py 


2.75 


32 - 152 


7 


Frascati 


S 


235 - 404 


90 


8 


Stanford 


Py 


172 - 436 


39 - 126 


19 


Bonn 


Py 


282 - 405 


74 - 98 


4 


Yale 


Py 


7 - 30 


48, 94 


20 


Yale 


Py 


7 - 13 


90 


3 


Tokyo 


Py 


352 - 697 


45 - 133 


27 


Kharkov 


S 


80 - 600 


75 - 150 


109 


Frascati 


S 


10 - 69 


90 


9 


Kharkov 


Py 


375 - 700 


43, 78, 90, 120 


40 


Kharkov 


Py 


550 - 1125 


90, 120 


30 


Yerevan 


S 


400 - 700 


45, 55 


5 


Tokyo 


T 


324 - 672 


72, 100, 130 


24 


Kharkov 




300 - 600 


75, 90, 120 


22,20,20 


Frascati 


S 


20, 29, 39, 61 


14 - 165 


41 


Bonn 


s 


233 - 818 


114, 135 


103 


Argonne 


Py 


6 - 14 


90 


6 


Kharkov 


E 


40, 50, 60, 70 


75, 90 


8 


Bonn 


T 


450, 550, 650 


25 - 155 


41 


Yerevan 


S 


395 - 795 


45 - 95 


30 


TRIUMPF 


An 

y 


180, 270 


32 - 144 


18 


TRIUMPF 


An 


370, 478 


45 - 155 


45 


Tomsk 


S 


50 - 100 


45, 60, 90 


13 


Wisconsin 


An 

y 


6, 13 


94 


2 


Novosibirsk 


T21 


33 - 125 


50 


4 


Yerevan 




300, 400, 500 


130 


3 


Kharkov 


Py 


200 - 367 


25 - 110 


30 


Yerevan 


Py ) Pxz 


306 - 436 


65, 75 


2x8 


Karlsruhe 


An 

y 


19 - 50 


62, 98, 131 


27 


Yerevan 


S 


284 - 999 


45, 60, 75, 95 


94 


BNL LEGS 


s 


100 - 314 


16 - 162 


112 


Novosibirsk 


T20 , 722 


49 - 505 


88 


2x9 


lUCF 


r' A A'^ 


183 


48 - 125 


3x6 


PSI 


An 

^y 


68 


69 - 144 


5 


Mainz 


S 


160 - 410 


35 - 155 


140 


Yerevan 


S 


787 - 1566 


90 


6 


JLab 


Pyi Cx', Cz' 


479 - 2411 


90 


10,9,8 



1971, 198f 
19 
20 
16 
16 
20 
20 



20q, 207, E0|, 20c, 21c ] 




22q , |22l|] 

i22|, m, |i7^, mi] 



T67|, 203, 204, 205] 



225| , |22q] 
22 
22 
22 



230, 231] 



23|, |23|] 



23|, m 

23l] 



237| , |238 
2p 

24 
24 

24^, |24| 



24|, |24|] 
2p] 
24 
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25 
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Figure 29. Four observables for d{'j,p)n at 300 MeV. Calculations are from Schwamb 



and Arenhovel (solid line) |175] and Kang, Erbs, Pfeil and RoUnik (dashed line) |177] 



The cross section data are from LEGS (triangles) |25C| and Mainz [p56| plus Bonn 



257 
254 



258, p59{ (circles). The S data are from LEGS (triangles) |250|, Mainz (circles) 



Kharkov (diamonds) ||210[, and Yerevan (stars) p47[. The T datum is from 



Tokyo [221|. The py data are from Stanford (triangles) | ]165[ , Bonn (circle) |166], 
Yerevan (star) |245| , and Kharkov (diamonds) ||22^ , |243| . See Table |ll| for related 
references. 



used three detector systems, two targets, and two different laser frequencies. Data were 
taken for = 100 - 315 MeV witli eiglit laboratory angles from 15 - 155°. Cross 
section statistical uncertainties range from a few percent to about 15%, with systematic 
uncertainties of 5%. One observation in the LEGS data is of pion contamination that 
may have been missed in earlier experiments, leading to increased cross sections. 

The MAMI Mainz experiment [p56|| used the Glasgow photon tagger along with the 
large solid angle detector DAPHNE to determine cross sections in the ranges E^ = 100 - 
800 MeV and 6cm = 30 - 160°. Data were binned by 20 MeV in energy and 10° in angle. 
Statistical uncertainties ranged from a few percent at lower energies to about 25% at 
the highest energies. Systematic errors were also energy dependent, ranging from a few 
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Figure 30. Four observables for d{^,p)n at about 450 MeV (410 MeV for S). 
Calculations are from Schwamb and Arenhovel (solid line) |175{ and Kang, Erbs, Pfeil 



and RoUnik (dashed line) |177|. The cross section data are from Mainz [ p56[ plus Bonn 
[ p57| , p58| . The S data are from Mainz (circles) p54| , Kharkov (diamonds) | plO{ , and 
Yerevan (stars) [247|. The T data are from Tokyo (square) [ p21[ and Bonn (circles) 

Tokyo (squares) |205| , Yerevan 



[231 . The py data are from Stanford (triangles) |16f: ^ ^ 

(star) m, and Kharkov (diamonds) l2l|, |22|. See Table |l| for related references. 



percent to several percent. The S asymmetry ||254|| was obtained by using a coherent 
bremsstrahlung radiator. 

Agreement between the Mainz and LEGS cross section results is generally better 
than 10%. The S asymmetries also agree well with each other and with earlier results 
from Yerevan ||247|| . Measurements from Kharkov [^D^, |207| , |209| , |21U|| generally agree, 
except for a tendency to be slightly smaller at many beam energies. 

The vector polarized target asymmetry T was measured at Tokyo [p20|, |221|| and 
at Bonn |230, 231|. The two measurements generally agree, with the data appearing to 



follow, very roughly, a — sin 6* dependence at each energy, as do the calculations of [p. 77 



The induced polarization py has been measured at a number of laboraties, with 
significant amounts of lower energy data from Stanford [|165|| , Tokyo ||167| , |203| , p04 
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205||, and Kharkov. There were several experiments at Kharkov, including initial 



measurements ||213| , |214| , |215|| , high-energy measurements ||216| , |217| , |218|| , simultaneous 



measurements of S, Ti, and Py ||222| , |223| , |170| , |224|| , and lower-energy measurements 
242l |243|| ■ The Kharkov data do not have a desirable level of consistency. The 



simultaneous measurements of E, Ti, and Py were taken as single-arm data, away from 
the photon endpoint, and may suffer from backgrounds. Polarizations below about 300 
MeV, including the intermediate energy neutron measurements ||234| , |235| , p36|| , appear 
to be well reproduced by theories. There are numerous data, particularly at 9cm = 
45°, 78°, 90°, and 120°, but there are few energies at which there are good angular 
distributions. The conclusion that polarizations are large, close to —1, and peak at 
about 500 MeV near 90° is beyond dispute. 



Figures ^ and ^ show the good agreement of the recent Mainz ||175|| and older 
Bonn [|1771] calculations with the cross section and E asymmetry. Theory seems to agree 
with da/dQ\\ better than da/dQ± ||249|| . The agreement is better at the lower energy, 
and the newer Mainz calculation is generally in better agreement. However, there is 
difficulty, particularly at the higher energy, with T and Py, both imaginary parts of the 
interference of amplitudes. The large induced polarizations above the A resonance have 
remained a puzzle for almost 30 years, and are still not fully explained by the newest 
theories. The Bonn calculation was in sufficiently good qualitative agreement with all 
observables {da/dfl, S, T, and py) for energies up several hundred MeV for the authors 
to consider this Py puzzle solved. However, detailed examination of Fig. ^ shows that 
neither the shape nor strength of the angular distribution is accurately reproduced; this 
will become clearer when we examine the energy dependence of Py at 6cm = 90° below. 



5.5.2. High energy photodisintegration Figure ^ shows the published high energy 
photo disintegration data, from experiments NE8 [p62| , p63|| and NE17 [|264|| at SLAC, 
and E89-012 |^ and E96-003 at JLab. These experiments determine cross sections for 
6cm ~ 36°, 52°, 69°, and 89° at energies from about 0.7 to 5.5 GeV; there are also some 
backward angle data up to 1.8 GeV from NE8. These data overlap well; the experiments, 
while all run by essentially the same collaboration, used three spectrometers in two 
experimental halls at two laboratories. There is also good overlap, variations of less 



than about 20%, with the highest energy Mainz tagged photon data ||256|| , and with 
older untagged data [gS^, |25|, |267| . 

Tagged photon measurements at low energies provide an accurate measure of beam 
flux, and along with the measured proton angle and energy, can determine a missing 
mass that allows background rejection. At high energies, smaller cross sections cannot 
be determined with the reduced flux of tagged photons. The Bremstrahlung endpoint 
technique was used for all of the SLAC and JLab measurements shown in Fig. [H^. 

In the endpoint technique, the measured proton momentum vector determines the 
incident photon energy and neutron kinematics, assuming the reaction is two body 
photodisintegration. Low momentum protons are cut from the analysis to prevent 
contamination from flnal states such as pnn^, while high momentum protons are cut 
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Figure 31. Photodisintegration cross section s^^da/dt versus incident lab photon 
energy. The calculations are from Kang, Erbs, Pfeil and RoUnik (solid line; 90° 



low energy data only) |177| , Lee (long dashed line; 90° low energy data only) |174[| , 
Raydushkin (dot-dashed Hne) RNA of Brodsky and Hiller (dotted) fTsl, quark 



qluon string model (short dashed) [190, 191|, and Frankfurt, Miller, Strikman and 



Sargsian (shaded region) |188, 189 



to eliminate the larger uncertainty in the photon flux close to the photon endpoint. 
Backgrounds are determined by radiator out and empty target measurements, and 
subtracted. Events in the region beyond the endpoint can be used to check the 
subtraction. The increase in time required for this subtraction makes it prohibitive 
for the highest energy measurements; for these the electrodisintegration background is 
calculated [ P6CI|| and subtracted. To determine the cross section, the incident photon 
flux is calculated using the method of Ref . ||261|| . Thick radiator corrections are typically 
about 15% for a radiator thickness of 6% of a radiation length. 
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The main feature of the cross section data above about 1 GeV is the s~^^ 
fall off of the cross sections da/dt {da/dQ) at = 90° and 69°, in agreement with 
the CCR and thus with perturbative QCD expectations. In contrast, the cross sections 
at the forward angles 36° and 52° fall off more slowly, with approximate scaling at 
lower energies until the onset of the s~^^ behavior at about 4 and 3 GeV beam energy, 
respectively. At each angle, the onset of the s~^^ behavior corresponds to ~ 1 GeV. 



The RNA and the Radyushkin estimates in Fig. |31| were normalized to the datum at 
89° and -E^ = 4 GeV, fixing their one free parameter. The RNA is then almost a factor 
of 2 too large at 36°, and also much too large at lower energies, requiring that the soft 
physics missing from the RNA interfere destructively with the leading terms. Suggested 



angular dependences of /(^cm) ||186|| would increase the RNA curve further at the forward 



angles, worsening agreeement with the data. In contrast, the Radyushkin estimate gives 
a somewhat better account of both the angular and energy dependence (even though it 
only goes asympotically as s~^^), confirming that phase space and nucleon form factors 
are all that is needed to account for much of the kinematic variation of the cross section. 
While the apparent onset of scaling at the forward angles suggests this agreement is 
starting to break down, we conclude the present data are insufficient to uniquely fix the 
asymptotic energy dependence of the cross section. 

The cross section data are also reasonably well reproduced by the model of 
Frankfurt, Miller, Strikman and Sargsian [FMSS] ||188| , |189|| and the quark gluon 



string (QGS) model ||190| , |191|| . The predictions of FMSS are uncertain because (a) 
the high energy NN scattering data has an uncertain energy dependence reflecting 
the experimental errors, and (b) the extrapolation of the NN data required for the 
predictions introduces further errors. These two uncertainties combine to give the 



jagged region shown in Fig. |31|. The QGS model describes the forward angle data 
up to 4 GeV reasonably well, even for values of —t exceeding the nominal limits 
of the model. The newer work predicts that the angular distributions will 



become increasingly symmetric at higher energies; older estimates [|190|] had predicted 
an increasing asymmetry . 

Two experiments currently have unpublished data for cross sections at photon 
energies up to 2.5 GeV. JLab Hall B E93-017 m used the CLAS with tagged photons to 



determine nearly complete angular distributions. The preliminary data agree well with 
earlier measurements, but are much more comprehensive than previous measurements 
in this energy range. Hall A E99-008 [p69|| has taken angular distributions at eight 



angles with 6'cm = 30° - 143°, at energies of 1.67, 1.95, and 2.50 GeV. Fig. ^ shows a 
sample angular distribution at approximately 1.6 GeV. 

5.5.3. Polarization observables in high energy photodisintegration There are only three 
sets of polarization data for deuteron photodisintegration at energies near and above 
1 GeV. The induced polarization py was measured at Kharkov ||216|, pl7| , |218|| , the E 



asymmetry was measured at Yerevan [p47| , p55|] , and Py and the polarization transfers 



and Cz' were measured at JLab |^. Data for the energy dependence of these four 




Figure 32. Deuteron photdisintegration angular distribution at ^ 1.6 GeV. Data 
are from SLAG experiments NE8 



262| and NE17 |264|, and JLab experiment 



E89-012 1^. Calculations are from Dieperink and Nagornyi [TtqI , Radyushkin |l8^ , 
and the RNA of Brodsky and Hiller |l86| . The Radyushkin curve has been normalized 
to the 4 GeV 90° datum, as in Fig. while the RNA curve has been reduced 60% 
from that normalization to better agree with this data. 



observables at a fixed 6cm = 90° are compared with theory in Figs. and 

The Kharkov measurements of Py at 6cm = 90° and 120° extend up to about 1.1 GeV 
(see Fig. |3^). These experiments were very difficuh. The small duty factor at Kharkov 
increases instantaneous background rates a factor of about 20,000 over those at JLab. 
These large backgrounds made it necessary to use multiple spark chambers to track 
particle trajectories. It was also difficult to calibrate the polarimeter. Calibrations of a 
polarimeter are best done by measuring its analyzing power using the known ep — >■ e'p 
elastic scattering reaction ^ ^ , but Kharkov had no polarized beam. The Kharkov 
measurements relied on a single elastic ep point to check false asymmetries in their 
polarimeter, and used analyzing powers from the literature. Finally, the polarimeter 
had a rear trigger scintillator; any inefficiencies in the scintillator would lead to false 
asymmetries. In contrast, the recent JLab experiment had little background, used 
ep — > e'p calibrations to determine false asymmetries and analyzing powers at each 
kinematic setting, and had no rear trigger scintillator. Given the clear disagreement 
of the Kharkov data with the recent JLab data, and noting that one of us (RG) is a 
spokesperson of the JLab experiment, we conclude that the highest energy set of Karkov 
data should not be trusted. 

The induced polarizations shown in Fig. ^ confirm our comments above concerning 
the Bonn calculation. The large negative polarization near 500 MeV is not reproduced, 
and the calculation is qualitatively incorrect at higher energies. The imaginary part of 
the amplitude appears to be a problem in these meson-baryon calculations; presumably 
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Figure 33. Induced polarization py for photodisintegration at 6'cm = 90°. Data 
sets from Stanford |165|, Tokyo |205|, and JLab are in the top panel. Data sets 
from Kharkov (squares) 



jl5[ , (stars) |218|, and (diamonds) |224| are in the lower 
panel. See Table [ll| for related references. The calculations are from Kang, Erbs, Pfcil 
and RoUnik (dash line) |177|, Sargsian (dash dot line) 270 and from Schwamb and 
Arenhovel (solid line) [175|. 



this arises from an inadequate treatment of resonances. 

Taken together, these recoil polarization data only weakly confirm the predictions 
of HHC, which predicts that Py and Cx' should approach zero as s — 00, and that 
(with additional assumptions ||178| , |185|| about relations between the helicity conserving 
amplitudes at 6'cm = 90°) S — 1 and C^i — >■ as s — > cxd. The highest energy 
polarization measurements of Py show that it is consistent with vanishing at energies 
above about 1 GeV, the same energy at which the s~^^ cross section scaling begins. [In 
the Radyushkin model Py should be zero because the amplitudes are all real if there is 



no gluon exchange ||271|| .] Similiarly, the polarization transfer observables Cx' and Czi 
both appear to peak near 1 GeV, and decrease at higher energies. However, Cx' does 
not appear to vanish sufficiently rapidly; the data might be inconsistent with HHC. So 
Py and Cz' (and perhaps Cx') seem to have close to the correct behavior, but S does 
not. The highest energy S asymmetry measurements from Yerevan, with data up to 
about 1.6 GeV, give the immediate impression that there is a minimum near 1.2 GeV 
and that above this energy the asymmetry is tending to increase towards 1, although 
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Figure 34. Polarization transfer for deuteron pliotodisintegration at Ocm = 90°. The 
data are from Ref. Tlie calculation is from Schwamb and Arenhovel |175 . HHC 
(with some additional assumptions for Cz') predicts that both of these amplitudes 
should vanish as s — > oo. 



the data are also statistically consistent with a constant value of about 0.3. In either 
case the trend is clearly not consistent with E — — 1. 

In models based on meson-baryon degrees of freedom, the data indicate that the 
combined effects of resonances plus final-state interactions are small. Calculations of Py 
in meson baryon theories ||171| , |172|, |173| , |176|, |177|1 indicate that Py at higher energies 
arises largely from resonance - background interference, with a small contribution from 
final-state interactions. Calculations generally indicate that the A resonance generates 
a large polarization, though only perhaps about 50% of the magnitude seen in the 
experimental data at = 90°. The Roper and Su have small effects, while the D13 
has a large effect ||176| , |177]| . The D33, included only by the Bonn group [ p.77|| , also 
generates a large polarization. (Of the 17 resonances included in the Bonn calculation 
177|| , only those mentioned above had large effects on py.) As discussed above, it is 



hard to imagine that a theoretically acceptable high energy model based on hadronic 
degrees of freedom will be constructed in the near future. Still, a modern relativistic 
calculation based on hadronic degrees of freedom would nonetheless be desirable. 

A calculation of Py has been done in the model of FMSS ||270|| . Since the helicity 
amplitudes in the nucleon-nucleon scattering for this center of mass energy range are 
not uniquely determined, some modelling was needed. The calculation showns that Py 
is generally very small, going from a negative value at low energies to positive values at 
several GeV beam energy, and is consistent with the trend of the experimental results. 
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Figure 35. Polarized plioton asymmetry S at ^cm = 90°. Tlie calculations are from 



Kang, Erbs, Pfeil and RoUnik (dot-dashed line) 177 and from Schwamb and Arenhovel 
(solid line) 



75 . The data are from Liu65 19 



Barbiellini67 IgOg, del BiancoSl pm , 
Gorbenko82 [|lO|, Adamian91 [|47|, and AdamianOO HI. See Table |ll| for related 
references. 



The polarization observable Cx is also expected to be small, and opposite in sign to Py, 
while Cz is expected to vanish. 



5.6. Future prospects 

The rapid falloff of photodisintegration cross sections with energy makes extension of the 
measurements difficult. As beam energy increases from 4 to 5, 6, and 7 GeV, the s~^^ 
dependence reduces cross sections by factors of 6, 30, and 115, respectively. Only a few 
experiments are possible without the proposed 12 GeV upgrade to JLab. As indicated 
in Table |l|, there are two approved experiments to continue the Hall A recoil polarization 
measurements to additional angles, and to beam energies near 3 GeV. Measurements of 
the S asymmetry are possible in JLab Hall B, over a range of angles and energies well 
over 1 GeV. 

In the longer term, the JLab 12 GeV upgrade offers additional possibilities. The 
luminosity increase planned for Hall B would allow precise polarization measurements 



to continue above 2 GeV. The proposed MAD spectrometer for Hall A ||272|| would 
give about a factor of 5 improvement in solid angle, and if it has the low backgrounds 
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characteristic of the current HRS spectrometers, cross section measurements are possible 
up to 7 GeV and polarization measurements are possible up to 4 GeV. 

5.7. Conclusions to Sec. 5 

Review of deuteron photodisintegration suggests the following: 

• A microscopic meson-baryon theory of deuteron photodisintegration must describe 
the NN interaction at high energies, including pion production and the 
contributions of hundreds of N* channels. It is unlikely that such a theory will 
be constructed in the foreseeable future. The data might indicate that the effect 
of many resonances is to increase the cross section and decrease the polarization 
obscrvablcs (by averaging over may phases) . This suggests that it might be possible 
to construct an effective theory based on hadronic degrees of freedom. 

• For p"^ greater than about 1 GcV^, cross sections appear to follow the constituent 
counting rules, but it is expected that an absolute pQCD calculation of the size of 
the cross section would give a result much too small. Similar observations may be 
made for other photoreactions, and it remains to be seen how this behavior arises, 
and if there is a general explanation for it. 

• The energy dependence of the photodisintegration cross sections has been shown to 
be potentially misleading indicator of the success of pQCD. Models with asymptotic 
behavior which differ from pQCD fit the data as well or better than pQCD. Further 
theoretical development and experimental tests of nonperturbative quark models 
would be desirable. 

6. Overall Conclusions 

The new high energy measurements of the deuteron form factors and the deuteron 
photodisintegration observables have motivated much theoretical work. In particular: 

• Conventional meson theory works well in cases where all of the active hadronic 
channels that can contribute to a process are included in the calculations. This has 
been done for the deuteron form factors (where only the NN channel is active) , but 
not for high energy deuteron photodisintegration where lOO's of N*N* channels are 
active. At high energy any successful meson theory must include relativistic effects. 

• New approaches, probably using quark degrees of freedom, are needed for high 
energy deuteron photodisintegration. While photodisintegration (as well as other 
reactions) seem to follow the scaling laws, theoretical estimates using pQCD give 
cross sections orders of magnitude too small. The data do not support hadronic 
helicity conservation. Thus scaling is no longer seen as sufficient evidence for the 
applicability of pQCD. 
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There is good evidence that, in this energy region, the deuteron is undergoing a 
transition from a region in which conventional hadronic degrees of freedom describe the 
physics to a region in which quark degrees of freedom are more appropriate. 
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